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ABSTRACT 


«& 


The  effects  of  pressure,  inlet  air  temperature,  and  fuel  type  on  the  soot 
threshold,  or  critical,  equivalence  ratio  are  presented.  Higher  pressures  yield 
lower  soot  thresholds,  while  no  dependence  on  fuel  type,  as  described  by 
either  the  fuel  hydrogen-to-carbon  ratio,  fuel  molecular  weight,  number  of 
carbon  atoms,  or  number  of  carbon-carbon  bonds,  is  observed.  Variations  in 
inlet  air  temperature  have  a  complex  effect;  however,  the  results  clearly  show 
that  the  experimentally  measured  flame  temperature  is  central  to  a 
description  of  the  incipient  soot  formation  process.  The  critical  equivalence 
ratio  dependence  on  pressure  and  temperature  is  shown  to  agree  with  a  two- 
step  semi-global  model  for  soot  precursor  evolution  for  pressures  from  0.1  to 
0.8  MPa,  and  measured  flame  temperatures  between  1600  and  2400K. 

The  effects  of  equivalence  ratio,  pressure,  and  fuel  chemistry  on  total 
non-luminous  flame  radiation  were  also  studied.  Radiant  intensity  was 
highest  for  an  equivalence  ratio  of  unity  and  increased  linearly  with  pressure 
from  0.4  to  0.8  MPa.  The  number  of  carbon-carbon  bonds  and  the  H/C  mass 
ratio  are  shown  to  influence  flame  radiation.  Empirical  non-luminous 
emissivities  increase  with  pressure  and  decrease  with  axial  location  and 
equivalence  ratio;  calculated  emissivities  show  no  trends. 


INTRODUCTION 


Soot  is  an  undesirable  product  of  combustion.  It  is  suspected  of  being 
harmful  to  human  health,  it  leaves  a  visible  exhaust  plume,  and  it  can 
overheat  combustor  walls  via  radiative  heat  transfer.  It  is  therefore 
advantageous  to  minimize  soot  formation  and  to  promote  oxidation  of  any 
soot  that  is  formed.  An  understanding  of  the  basic  relationship  between 
operating  conditions  and  soot  formation  is  necessary  to  accomplish  this.  In 
particular,  this  study  undertook  as  a  first  major  objective  the  determination 
of  the  dependence  of  soot  formation  on  equivalence  ratio,  pressure, 
temperature,  and  fuel  being  burnt. 

Flame  radiation  is  a  second  important  phenomenon  associated  with  the 
combustion  process.  In  some  applications,  for  example  in  a  power  plant 
furnace,  high  radiative  wall  heat  fluxes  are  desirable  [Wark,  1977].  In  other 
systems  radiative  heat  loads  can  overheat  vital  components.  Excessive 
radiative  heat  flux  is  particularly  critical  in  gas  turbine  combustion  chambers 
because  film  cooling  of  the  combustor  liner  provides  an  effective  barrier 


against  convective  heating  but  not  against  radiative  transfer. 

The  gas  turbine  combustor  has  been  used  as  a  test  bed  in  previous  flame 
radiation  studies  at  high  pressures;  however,  those  studies  have  all  used 
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multicomponent,  heavy-hydrocarbon  liquid  fuels.  Consequently,  they  were 
unable  to  isolate  the  effects  of  fuel  structure  on  flame  radiation  at  high 
pressures.  This  study  therefore  undertook  as  a  second  major  objective  an 
investigation  into  flame  radiation  from  gas-fueled  flames. 


The  overall  goal  of  this  investigation  is  to  extend  the  fundamental 
knowledge  of  soot  formation  and  flame  radiation  to  combustion  at  elevated 
pressures.  The  principal  application  of  this  information  lies  in  the  design  of 
high-pressure  combustors  for  minimum  soot  emission  and  optimum  radiative 
heat  transfer  characteristics. 


Soot  Formation 


The  detailed  chemical  mechanism  for  soot  formation  is  not  well 
understood  but  can  be  modeled  in  a  semi-global  manner  as  a  two-stage 
process  [Harris  et  al,  1886].  The  first  step  occurs  when  fuel  molecules  enter 
the  pyrolysis  zone  of  a  flame  and  are  converted  to  a  variety  of  compounds, 
including  soot  precursors.  The  soot  precursors  are  subsequently  converted  to 
either  soot  or  oxides  of  carbon,  depending  upon  the  nature  of  the  flame. 
When  the  formation  of  soot  precursors  occurs  at  the  same  rate  as  their 
oxidation,  the  precursor  concentration  is  at  steady  state,  and  the  flame  is  at 
the  soot  threshold.  The  equivalence  ratio  at  this  point  is  termed  the  critical 
equivalence  ratio,  <f>c. 


The  soot  threshold  itself  is  an  important  experimental  quantity  for  several 


reasons.  Since  it  is  the  fuel/air  mixture  at  which  the  rates  of  soot  precursor 
formation  and  destruction  balance  each  other,  its  measurement  can  provide 
information  to  determine  the  extent  of  applicability  of  the  two-stage  model 
described  above.  The  measurement  is  very  simple  and  repeatable,  depending 
only  upon  the  visual  acuity  of  the  observer.  Finally,  the  soot  threshold 
provides  a  quantitative  measure  of  the  onset  of  flame  luminosity,  obviously 
important  when  analyzing  radiation  characteristics. 

Several  important  questions  about  soot  thresholds  remain  to  be  answered. 

—  Is  the  two-step  global  model  valid  for  combustion  pressures  above  one 
atmosphere? 

—  Which  flame  temperature  best  describes  the  sooting  condition:  adiabatic 
or  measured?  And  which  measured  temperature? 

—  Will  an  effective  equivalence  ratio,  based  on  products  of  carbon  monoxide 
and  water  vapor,  describe  results  more  accurately  than  the  standard 
equivalence  ratio,  based  on  products  of  carbon  dioxide  and  water  vapor? 

—  How  important  are  fuel  properties,  such  as  molecular  weight,  hydrogen- 
to-carbon  ratio,  and  number  of  carbon-carbon  bonds? 

This  study  answers  these  questions  by  investigating  soot  thresholds  at  high 
pressures. 

Since  the  experiments  in  this  work  were  run  at  various  temperatures  and 
pressures  and  since  the  semi-global  model  requires  the  determination  of  OH 


concentrations,  a  knowledge  of  equilibrium  OH  behavior  with  respect  to  these 
variables  is  required.  For  instance,  an  increase  in  temperature  brings  about 
an  exponential  increase  in  the  OH  concentration.  Higher  pressures,  on  the 
other  hand,  lead  to  decreases  in  the  OH  concentration.  Thus,  from  these 
considerations  we  should  expect  high  flame  temperatures  and/or  low 
pressures  to  preclude  sooting.  The  pyrolysis  of  the  fuel,  however,  is  also 
temperature  dependent,  higher  temperatures  giving  higher  concentrations  of 
soot  precursors.  This  effect,  then,  counteracts  the  higher  OH  concentrations 
brought  about  by  higher  temperatures.  A  significant  contribution  to  the 
understanding  of  soot  formation  could  be  made  if  the  two-stage  model 
described  above  could  be  shown  to  describe  the  critical  equivalence  ratio 
behavior  found  in  this  investigation. 


Flame  Radiation 


The  thermal  radiation  spectrum  encompasses  wavelengths  between  0.4 
and  1000  /im  and  includes  band  as  well  as  continuous  radiation.  Band 
radiation  is  emitted  by  gaseous  combustion  products,  principally  carbon 
dioxide  and  water  vapor.  The  strongest  carbon  dioxide  bands  are  centered  at 
2.7  /im  and  4.3  /im,  while  the  water  bands  are  located  at  1.9  /im,  2.8  /im, 
and  6.7  /im.  Continuous  radiation  emitted  by  hot  soot  particles  extends  over 
all  thermal  wavelengths.  However,  even  for  a  cool  flame  of  1600  K,  95%  of 
the  energy  is  emitted  at  wavelengths  shorter  than  8.0  /im.  For  a  flame  at 


2000  K  the  95%  cut-off  is  6.3  urn.  The  significant  conclusion,  then,  is  that  by 


measuring  radiative  energy  at  wavelengths  between  0.7  /um  and  8.0  £zm  one 
observes  essentially  all  flame  radiation  pertinent  to  this  study. 

Radiation  is  important  in  this  study  because  it  is  a  major  contributor  to 
the  heat  flux  at  the  combustor  wall.  Radiation  is  a  complex  phenomenon 
and  is  only  partially  understood;  nevertheless,  the  following  parameters  are 
known  to  affect  the  wall  heat  flux  [Lefebvre,  1983].  Combustor  pressure  can 
alter  soot  levels  and  emission  band  widths,  thereby  altering  radiation  fluxes. 
Flame  temperature  strongly  influences  flame  radiation  and  is  controlled  by 
equivalence  ratio  and  inlet  air  temperature.  Equivalence  ratio  and 
differences  in  fuel  chemical  structure  can  lead  to  changes  in  the 
concentrations  of  C02  and  H20,  upon  which  flame  radiation  also  depends. 
For  these  reasons  this  study  focuses  on  measurement  of  thermal  radiation 
emitted  by  flames  at  various  values  of  pressure,  inlet  air  temperature, 
equivalence  ratio,  and  fuel  type. 

The  studies  of  flame  radiation  at  high  pressures  are  numerous  [Moses  and 
Naegeli,  1979;  Naegeli  et  al.,  1983;  Claus,  1981;  and  Humenik  et  al.,  1983], 
but  the  results  have  been  influenced  by  the  interdependences  of  the  variables 
involved.  The  goal  of  the  radiation  portion  of  this  study  is  to  isolate  the 
effects  of  the  various  parameters  studied.  A  particularly  important  aspect, 
then,  is  that  the  fuels  be  burned  in  the  gaseous  state  thereby  eliminating  any 
complicating  effects  of  fuel  spray  distribution  and  drop  size,  a  problem 
common  in  liquid-fueled  combustors.  Furthermore, -o»ly  pure  hydrocarbon 


fuels  should  be  burned.  The  effect  of  fuel  chemical  structure  is  thus  more 
evident  than  in  studies  using  multi-component  fuels. 

The  flame  radiation  intensity,  If,  is  the  radiant  power  emitted  per  unit 
area  projected  normal  to  the  direction  of  passage,  per  unit  solid  angle.  It  is 
the  radiometric  quantity  reported  in  this  study.  Because  the  field  of  view  of 
the  radiometer  is  centered  on  a  diameter  of  the  flame,  the  emitting  area  and 
the  projected  areas  are  the  same.  This  is  equivalent  to  stating  that  the 
numerical  values  of  the  normal  emitted  flux  and  the  intensity  normal  to  the 
flame  are  interchangeable  when  using  this  radiometric  alignment.  The 
significant  point,  however,  is  that  these  quantities  were  measured  at  the 
boundary  of  the  cylindrical  flame.  As  such,  no  consideration  was  given  to 
attenuation  by  absorption  or  scattering  or  to  inward  scattering.  The 
reported  intensities  are  simply  the  bulk  quantities  emitted  normal  to  the 
flame  (i.e.  along  a  diametral  path)  as  measured  from  outside  the  boundary  of 
the  hot  flame  gases. 

The  radiation  intensity  reported  here  is  of  value  to  the  combustion 
engineer  as  a  design  guideline.  Many  practical  combustors  burn  gaseous 
fuels,  most  of  them  among  the  fuels  studied  here.  An  estimate  of  flame 
temperature  and  a  knowledge  of  combustor  geometry  are  all  the  information 
required  to  apply  the  results  given  in  this  work  to  an  actual  system. 


LITERATURE  REVIEW 


Critical  Equivalence  Ratio 

Recent  studies  of  soot  thresholds  in  premixed  flames  have  used  the 
minimum  equivalence  ratio  at  which  soot  was  detected  as  the  principal 
indicator  of  a  fuel’s  tendency  to  form  soot.  This  so-called  critical  equivalence 
ratio,  4>c,  is  most  often  based  on  complete  combustion  to  carbon  dioxide  and 
water  vapor.  Other  end  conditions  have  been  considered,  however. 
Takahashi  and  Glassman  [1984]  suggested  an  effective  equivalence  ratio  that 
assumed  combustion  products  of  carbon  monoxide  and  water  vapor.  Olson 
and  Pickens  [1984]  and  Millikan  [1962]  discussed  the  use  of  the  carbon-to- 
oxygen  ratio  in  the  unburnt  mixture,  based  on  combustion  products  of  carbon 
dioxide  and  hydrogen.  However,  Olson  and  Pickens  state  that  none  of  these 
effective  equivalence  ratios  is  superior  to  the  others,  probably  because  the 
combustion  products  vary  so  greatly  for  different  equivalence  ratios.  Because 
the  definition  of  equivalence  ratio  based  on  carbon  dioxide  and  water  vapor 
is  the  most  widely  known,  it  will  be  used  here.  (Note  that  a  high  value  of  <f)c 
indicates  a  low  tendency  to  form  soot.) 
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Effect  of  Temperature 


Flame  temperature  has  been  widely  studied  as  a  parameter  having  an 
effect  on  <j)c.  Glassman  and  Yaccarino  [1981]  found  a  distinct  relationship 

between  4>c  and  — .  Millikan  found  a  linear  relationship  between  his  effective 

critical  equivalence  ratio  (based  on  the  carbon-to-oxygen  atom  ratio  in  the 
unburnt  mixture)  and  measured  flame  temperature  for  ethylene/air  flames. 
He  also  showed  a  linear  relationship  between  4>c  and  inverse  temperature.  He 
suggested  that  ln[OH]  was  a  linear  function  of  equivalence  ratio  and  that 
higher  temperatures  led  to  higher  [OH],  which  in  turn  led  to  more  rapid 
oxidation  of  soot  precursors  within  the  pyrolysis  zone. 

Harris  et  al.  extended  this  theory,  arriving  at  a  slightly  different 

<f>c 

expression.  They  found  a  linear  relationship  between  In- - -  and  inverse 

[OH] 

temperature  that  held  for  five  different  aliphatic  fuels  and  postulated  that 
the  formation  of  soot  precursors  was  first  order  in  fuel  concentration  while 
the  oxidation  of  the  precursors  was  first  order  in  [OH].  They  reasoned  that  a 
two-step  model— precursor  formation  followed  by  oxidation— would 
satisfactorily  predict  soot  thresholds.  They  also  postulated  that  at  the 
critical  equivalence  ratio  the  rate  of  formation  of  soot  precursors  exactly 
balanced  the  rate  of  destruction  of  precursors,  or 


[P]  =  kf[F]  -  kd[OH][P]  =  0 


where 


[P]  =  time  rate  of  change  of  the  precursor  concentration 

[F]  =  fuel  concentration 

[OH]  =  hydroxyl  radical  concentration 

[P]  =  precursor  concentration 

kf  =  rate  coefficient  for  precursor  formation 

kd  =  rate  coefficient  for  precursor  destruction 


Over  a  given  narrow  range  of  pressures  and  flame  temperatures  the  fuel 
concentration  is  a  linear  function  of  equivalence  ratio,  so 


[F]  -  c -4>  (2) 

With  the  reaction  rate  coefficients  taking  an  Arrhenius  form,  Eqn  (1)  can  be 


reduced  to 


E»,f  ~  E»,d 


where 


Ea  =  activation  energy  (f  =  formation;  d  =  destruction) 
Cfuej  =  a  constant  specific  to  a  given  fuel 


An  unexpected  result  came  from  the  work  of  Olson  and  Madronich  [1985], 
who  measured  temperatures  and  soot  threshold  equivalence  ratios  in 


aromatic-fueled  flames  as  functions  of  the  ratio  of  oxygen  to  nitrogen  in  the 


unburnt  mixture.  While  they  concur  with  Takahashi  and  Glassman  that  4>c 
increases  with  adiabatic  flame  temperature,  they  found  tnat  measured  flame 
temperatures  were  unrelated  to  soot  threshold  equivalence  ratios.  In  a 
similar  study  Wright  [1969]  adjusted  the  flame  temperature  within  a  well 
stirred  reactor  by  preheating  the  fuel/air  mixture.  Burning  olefins  and 
aromatics  with  air  preheated  to  as  high  as  650  K,  he  demonstrated  that  an 
increase  in  pre-flame  temperature  gave  a  higher  soot  threshold  equivalence 
ratio.  A  similar  rise  in  soot  threshold  followed  the  introduction  of  diluents  of 
lower  heat  capacity,  an  effect  which  he  also  attributed  to  an  increase  in 
flame  temperature.  His  experimental  method  had  an  important  difference, 
however,  from  all  others  cited  here;  his  criterion  for  critical  soot  limit 
required  his  mixtures  to  be  rich  enough  to  leave  soot  in  the  end  products  of 
the  flame.  All  other  soot  threshold  measurements  have  been  made  at  the 
point  of  incipient  soot  formation  where  even  though  soot  may  be  forming  at 
one  point  within  the  flame  it  can  be  burnt  off  in  downstream  locations. 
Thus,  Wright's  critical  equivalence  ratio  values  are  expected  to  be  higher 
than  those  of  the  other  studies  discussed  here.  Wright’s  pressure  and 
temperature  dependencies  may  also  differ  from  those  observed  at  the  critical 
soot  point. 

Most  of  these  studies,  then,  have  concluded  that  a  higher  flame 
temperature,  either  measured  or  adiabatic,  leads  to  a-hifber  4>c;  however,  the 


precise  relationship  between  flame  temperature  and  soot  threshold  needs 
refinement.  In  addition,  the  dependence  of  4>c  on  pressure  is  not  currently 
known. 

Effect  of  Pressure 

Pressure  effects  on  4>c  have  been  investigated  less  extensively  than 
temperature  effects.  Calcote  and  Miller  [1978],  burning  acetylene,  benzene, 
and  cumene,  found  that  soot  threshold  equivalence  ratios  were  independent 
of  pressure  between  9  and  100  kPa.  MacFarlane  et  al.  [1975]  also  determined 
the  soot  threshold  to  be  independent  of  pressure  over  the  range  0.5  to  2.0 
MPa  for  gaseous  pentanes  and  hexanes.  In  contrast,  Fenimore  et  al.  [1957] 
observed  a  decrease  in  (f>c  as  the  pressure  rose. 

The  resuks  of  Calcote  and  Miller  and  MacFarlane  et  al.  [1975]  are 
surprising  since  the  concentration  of  OH  radicals  is  known  to  be  a  function  of 
pressure.  Since  OH  radicals  are  believed  to  play  an  important  role  in  the 
oxidation  of  soot  precursors,  one  would  expect  that  as  pressure  changes,  [OH] 
would  change.  Based  on  the  model  of  Harris  et  al.,  the  soot  threshold  would 
also  change.  No  experimental  data  exist,  however,  to  either  verify  that  <f>c 
does  depend  on  pressure,  or  to  quantify  that  dependence.  For  this  reason  the 
effect  of  pressure  was  investigated  in  the  present  study. 


Effect  of  Fuel  Type 


Fuel  effects  on  4>c  have  been  the  subject  of  many  investigations.  Several 
prior  studies  have  been  performed  with  liquid  fuels,  but  in  many  cases  the 
results  are  pertinent  here.  Naegeli  et  al.  [1983],  for  example,  showed  that  the 
hydrogen-to-carbon  ratio  of  a  liquid  fuel  was  a  good  indicator  of  a  fuel’s 
tendency  to  soot.  An  important  corollary  was  that  the  internal  structure  of 
the  fuel  molecule  was  unimportant.  These  authors  also  pointed  out  that 
improved  aeration  of  a  flame  reduces  the  dependence  of  soot  on  flame 
temperature,  implying  that  premixed  flames  will  be  less  sensitive  to 
temperature  than  diffusion  flames.  Calcote  and  Manos  [1983],  in  an 
exhaustive  review  of  atmospheric  pressure  flames  using  a  variety  of  fuels  from 
several  hydrocarbon  families,  listed  several  important  trends  for  premixed 
flames.  First,  the  result  of  changing  the  carbon-to-hydrogen  ratio  from  fuel 
to  fuel  appeared  to  be  ambiguous.  All  alkenes,  for  example,  have  a  C/H 
value  of  l/2;  however,  the  soot  thresholds  for  these  fuels  varied  widely. 
Furthermore,  acetylene  (with  its  high  C/H  value  of  1)  had  the  highest  soot 
threshold  of  any  fuel  examined.  Second,  soot  thresholds  fell  as  the  molecular 
weight  rose  within  a  hydrocarbon  family  (e.g.  alkane  family).  This  was 
especially  apparent  for  hydrocarbons  of  fewer  than  six  carbon  atoms. 
Finally,  the  effect  of  molecular  structure  showed  the  following  trend  for  soot: 


highest  thresholds  (least  likely  to  soot) 
*  alkynes  (acetylenes) 

alkenes  (olefins) 
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iso-alkanes 
n-alkanes  (paraffins) 
lowest  thresholds  (most  likely  to  soot). 

Takahashi  and  Glassman  defined  a  new  parameter  for  describing  critical 
equivalence  ratios  in  premixed  flames,  namely  the  total  number  of  carbon- 
carbon  bonds  within  a  molecule.  They  showed  that  for  a  wide  variety  of 
aliphatic,  alicyclic,  and  aromatic  fuels,  all  burned  at  the  same  calculated 
adiabatic  flame  temperature,  the  critical  equivalence  ratio  decreased  as  the 
number  of  carbon-carbon  bonds  increased.  In  an  interesting  discussion  on 
the  effects  of  flame  temperature  and  fuel  type  Calcote  and  Olson  [1982] 
stated  that  "although  increasing  the  flame  temperature  of  an  individual  fuel 
may  decrease  its  tendency  to  soot,  a  series  of  different  fuels  with  increasing 
flame  temperatures  at  their  soot  points  generally  increase  in  tendency  to  soot 
(ethylene  and  acetylene  are  exceptions)."  The  important  point  here  is  that 
fuel  structure  may  be  the  dominant  factor  in  determining  the  soot  threshold. 

To  facilitate  comparison  of  data  from  all  soot  threshold  investigations 
Calcote  and  Manos  have  devised  a  Threshold  Soot  Index  (TSI)  for  premixed 
flames: 
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■where  4>c  is  the  soot  threshold  equivalence  ratio  and  a  and  b  are  constants 
specific  to  a  given  experimental  system.  The  goal  of  the  TSI  is  to  remove 
variations  from  rig  construction,  operational  procedures,  and  systematic 
experimental  variations,  thereby  bringing  together  the  results  of  many 
different  authors.  The  constants,  a  and  b,  are  computed  by  taking  (fic  data 
for  two  fuels  of  very  different  molecular  weights  and  fitting  them  to  the  now 
widely  published  TSI  values.  A  high  value  of  TSI  reflects  a  greater  tendency 
to  soot  (i.e.  low  <f>c),  while  a  low  value  of  TSI  indicates  a  lesser  tendency  to 
soot  (i.e.  high  4>c). 

The  present  study  has  two  experimental  characteristics  that  could  cause 
the  results  to  deviate  from  those  of  other  investigators.  First,  the  soot 
threshold  is  determined  visually,  subjecting  the  results  to  the  sensitivity  of 
the  human  eye.  Second,  the  flame  is  subjected  to  significant  cooling, 
especially  from  the  water  jacket  surrounding  the  combustor.  Neither  of  these 
characteristics  is  expected  to  diminish  the  precision  of  the  study,  but  the 
results  were  different  from  those  of  flat  flame  burners,  which  in  some  cases 
approach  adiabatic  performance. 


Flame  Radiation 


Little  radiant  heat  transfer  data  exist  for  high-pressure,  gaseous-fuel 
flames;  however,  many  studies  have  been  conducted  on  flame  radiation  from 
high-pressure  combustors  burning  liquid  fuels.  The  expression*  derived  from 
liquid-fuel  data  are  included  here  in  hopes  of  providing  an  insight  into  the 
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results  of  the  present  work.  The  complex  interactions  of  fuel  atomization 
and  combustor  geometry  with  flame  radiation  should  be  borne  in  mind, 


however,  in  any  comparison  between  these  dissimilar  experiments. 

Effect  of  Temperature 

An  analytical  method  developed  by  Kunitomo  and  Kodama  [1974] 
predicts  that  the  contributions  to  flame  emissivity  from  soot  and  non- 
luminous  gasses  are  completely  independent  of  flame  temperature  for  any 
rich  fuel/air  mixture  at  pressures  of  0.3  and  0.5  MPa.  No  experimental 
evidence  has  been  found  to  support  this  result. 

Effect  of  Pressure 

Several  investigators  have  examined  the  effect  of  pressure  on  flame 
radiation  in  practical  gas  turbine  combustors  [MacFarlane  et  al.,  1975; 
Marsland  et  al.,  1975;  Humenik  et  al.,  1983].  All  concur  that  higher  pressures 


yield  considerably  higher  radiative  fluxes  at  the  combustor  walls  with  a 
diminishing  effect  at  very  high  pressures  (1.2  to  1.5  MPa).  The  work  of 
MacFarlane  et  al.  [1975]  is  especially  relevant  since  he  measured  flame 
radiation  and  emissivity  for  premixed  methane/air  flames.  Using  the 
Schmidt  technique  he  found  that  the  radiative  flux  also  rose  with  pressure  for 
rich  fuel/air  mixtures.  In  addition,  he  found  that  an  increase  in  pressure 
brought  about  a  substantial  increase  in  emissivity,  especially  for  equivalence 
ratios  above  2.0.  While  this  information  may  be  useful,  discretion  is 
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EXPERIMENTAL  SYSTEM 


This  chapter  describes  the  hardware  and  instrumentation  included  in  the 
gaseous-fueled  combustion  rig.  It  begins  with  the  air  supply  system,  follows 
with  the  fuel  delivery  system,  and  then  discusses  the  combustor  as  a  whole. 
Finally,  the  three  combustor  subsystems  (igniter,  flameholder/fuel  injector, 
and  flame  positioning  mechanism)  are  discussed. 

Air  Supply  System 

The  laboratory  air  supply  system  consists  of  two  compressors  connected 
to  four  reservoir  tanks  (see  Figure  l).  Tank  refilling  can  proceed  during  an 
experimental  run,  if  desired.  A  mechanical  description  of  the  compressors, 
the  air  delivery  system,  and  the  air  regulators  is  presented  below.  A 
discussion  of  the  uses  of  the  four  air  subsystems,  classified  by  their  maximum 
working  pressure,  follows.  Finally,  a  description  of  the  air  flow  profile 
through  the  combustor,  and  its  measurement,  is  presented. 
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Mechanical  Components 

The  air  is  raised  to  its  dry,  high-pressure  state  using  a  set  of  two  piston- 
type  compressors  and  two  drying  systems.  The  first  set  of  compressors  boosts 
ambient  air  to  0.345  MPa  and  533  K.  The  air  then  passes  through  a 
separator  to  remove  liquid  water  and  oil.  Next,  the  air's  dew  point  is 
reduced  to  236  K  in  a  chamber  charged  with  activated  aluminum.  Dust  is 
then  removed  from  the  air  with  a  paper  filter  and  the  air  is  pumped  to  15.2 
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MPa.  Finally,  an  activated  carbon  filter  removes  hydrocarbons  from  the  air 
before  it  is  sent  to  the  reservoir  tanks. 
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Combustion  Air  Pre-T reatment 

Air  entering  the  combustor  was  heated  by  means  of  a  Wellman  model 
W15000  immersion  heater,  comprised  of  three  resistively  heated  cal  rods 
powered  by  a  240  V,  three-phase  SCR  controller.  At  the  maximum  current  of 
60  A,  14.4  kW  of  power  was  available.  A  4-20  mA  controller,  located  in  the 
control  room,  regulated  the  current  delivered  to  the  heater.  Its  reference 
signal  came  from  a  thermocouple  located  at  the  combustor  inlet. 

Incoming  air  supply  flow  conditions  were  carefully  maintained  during' 
these  experiments.  The  combustor  entrance  incorporated  two  bends,  • 
installed  to  fit  the  rig  within  the  test  cell  confines.  These  turns  distorted  the 
air  flow  field  causing  asymmetry  of  the  velocity  profile  and  imparting  swirl; 
however,  a  set  of  flow  straighteners  and  screens,  located  directly  upstream  of 
the  combustor,  re-established  uniform  flow.  The  length-to-diameter  ratio  of 
the  flow  straightener  was  44  to  1,  enough  to  provide  au  even  velocity  profile 
and  to  minimize  swirl  [Daily  and  Harleman,  1966].  A  set  of  60  x  60  mesh, 
0.19  mm  diameter  wire  screens  were  silver  soldered  upstream  of  the  flow 
straightener  in  accordance  with  the  work  of  Schubauer  et  al  [1948].  For  an 
inlet  temperature  of  600  K,  a  nominal  velocity  of  6  m/s  through  the 
combustor  gave  ReD  =  4300,  far  above  the  critical  ReD  of  2300  for  fully 
turbulent  flow.  Because  the  screens  were  placed  in  the  flow  and  because  the 
Reynolds  number  was  well  above  the  critical  value  for  turbulence,  all  flames 
were  considered  to  be  fully  turbulent. 


id 


Axial  air  velocity  at  several  mass  flow  rates  was  measured  with  a 
traversing  pitot-static  probe  connected  to  an  inclined  manometer.  The 
velocity  at  any  radial  position,  within  a  given  plane  normal  to  the  flow,  was 
found  to  lie  within  four  percent  of  the  average  value  across  this  plane.  See 
Figure  2.  All  velocity  measurements  were  made  with  the  fuel 
injector/flameholder  apparatus  removed  from  the  combustor. 

Fuel  Delivery  System 

The  fuel  system  was  designed  to  deliver  gaseous  fuels  to  the  combustor 
and  to  control  and  measure  their  flow  rates  (see  Figure  3).  Most  fuels  had  to 
be  heated  or  pressurized,  as  described  below,  to  propel  them  into  the 
combustion  chamber.  This  was  accomplished  by  either  warming  the  storage  . 
bottles  to  boost  vapor  pressures  or  driving  the  fuels  in  liquid  form  using 
nitrogen  until  they  could  be  vaporized  by  an  immersion  heater.  Three 
different  modes  of  fuel  storage  were  used:  (l)  Propane  was  stored  as  a  two- 
phase  mixture  in  a  large  bottle  outside  the  laboratory  building.  (2)  Butane, 
butylene,  cyclopropane,  and  propylene  were  stored  in  the  fuel  room  as  two- 
phase  mixtures  within  small  bottles.  (3)  Methane,  ethane,  and  ethylene  were 
stored  as  gases  in  standard  cylinders. 

Pressurization 

All  fuels  except  methane,  ethane,  and  ethylene  had  to  be  pressurized  to 
force  them  into  the  high-pressure  combustor. 

The  volume  above  the  liquid  in  the  propane  storage  tank  was  pressurized 
with  facility  nitrogen.  A  bottom-fed  line  tapped  the  pressurized  liquid 
propane  for  delivery  to  the  lab.  The  propane  bottle  was  fitted  with  a  1.72 
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MPa  relief  valve  to  prevent  overpressurization.  See  Figure  4. 

Bottles  containing  butane,  butylene,  cyclopropane  and  propylene  were 
placed  in  a  water  bath  (T  =  70  C)  to  boost  their  vapor  pressures.  This 


technique  allowed  cyclopropane  to  be  injected  at  pressures  up  to  0.4  MPa 
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and  propylene  up  to  0.6  MPa.  Butane  and  butylene  could  be  injected  at 
pressures  up  to  0.3  MPa. 

Fuel  Pre-Heating 


The  fuel  being  tested  was  heated  to  the  vapor  state  before  being  routed 
to  the  fuel  injector.  The  fuel  heater  design  was  similar  to  that  of  the  air 
heater,  but  was  powered  by  a  single-phase  SCR  controller  and  oriented 
vertically.  Maximum  heater  power  was  1.2  kW  and,  as  with  the  air  system, ' 
the  heating  rate  was  controlled  using  a  4-20  mA  temperature  monitor  that 


received  data  from  a  thermocouple  mounted  in  the  fuel  delivery  tube.  All 
fuels  were  heated  to  375  K,  a  temperature  sufficient  tc  vaporize  any  of  the 
fuels  used  here  without  pyrolysis  occurring  but  safely  below  the  melting 
temperature  of  the  buna  O-rings  present  in  the  fuel  supply  lines. 

Since  the  fuel  did  not  pass  immediately  from  the  immersion  heater  to  the 
injector,  care  had  to  taken  to  keep  it  in  the  vapor  state.  This  was 
accomplished  by  means  of  an  electrical  pipe  heater  (3  W/m)  wrapped  around 


the  delivery  line  from  the  heater  to  the  injector.  Fiberglass  insulation  was 
then  wrapped  around  the  heating  strip/fuel  delivery  line  to  optimize  heat 
transfer. 
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Control . 


The  fuel  regulating  valve  was  located  less  than  two  meters  from  the  fuel 
injection  tube  to  minimize  the  lag  between  valve  actuation  and  a  change  in 
combustor  performance.  Precise  control  of  fuel  flow  rates  was  facilitated  by 
two  component  modifications.  First,  the  regulating  valve,  originally  oversized 
for  the  fuel  flow  rates  intended  for  this  work  (1.0-20.0  g/s),  was  recut  in  the 
pintle/seat  chamber  to  give  it  a  very  gentle  taper.  Second,  remote  control, 
made  possible  by  using  a  Grove  hand  loader  delivering  69  to  86  kPa  pressure 
to  the  valve  diaphragm,  was  made  more  precise  by  replacing  the  hand  loader 
spring  that  came  from  the  factory  with  another  of  lower  spring  constant. 
Finally,  fuel  was  routed  through  flexible  hose  between  the  on/off  valve  and. 


the  injector  tube,  thus  permitting 


the  section  on  Flame  Positioning). 


Fuel  Purge 

Nitrogen  was  forced  through  the  fuel  heater  and  into  the  combustor  to 
purge  those  lines  of  any  residual  combustibles.  Actuation  came  from  a  28  V 
DC  on/off  valve  placed  just  upstream  of  the  fuel  heater,  protected  from  fuel 
back  flow  by  a  check  valve. 

Nitrogen  was  also  used  during  the  warm-up  of  the  fuel  heater.  Because 
the  controlling  thermocouple  for  the  heater  lay  just  downstream  of  its  outlet, 
gas  had  to  be  flowing  through  the  line  so  that  the  temperature  could  be  read 
properly.  Using  fuel  would  have  been  wasteful  and  dangerous,  so  a  small 
flow  of  nitrogen  was  used  instead. 


Combustion  Chamber 


The  combustor  consisted  of  a  5  cm  OD,  3.8  cm  ID,  XX-347  stainless  steel 
pipe,  cooled  by  two  water  systems:  (1)  a  spiraling,  low-pressure,  high-flow¬ 
rate,  external  jacket  where  the  water  was  confined  to  an  annulus  between 
the  external  wall  of  the  combustor  pipe  and  the  inner  wall  of  the  jacket  and 
(2)  direct  injection  of  water  downstream  of  the  windows  to  quench  the  flame, 
cool  most  of  the  length  of  the  suction  thermometer  and  keep  the  back¬ 
pressure  valve  cool  (see  Figure  5).  The  high  pressure  water  for  the  direct 
injection  cooling  system  was  provided  by  a  centrifugal  pump  driven  by  an 
electric  motor.  The  pump  was  plumbed  in  parallel  with  the  low-pressure 
water  supply  so  that  a  hand  loaded  valve,  when  closed,  forced  water  from 
the  low-pressure  branch  to  the  pump. 

Two  diametrically  opposed  windows  were  installed  in  the  combustor  to 
permit  measurement  of  radiative  flux  and  observation  of  soot  thresholds. 
Three  items  were  considered  when  designing  these  windows.  (1)  The  windows 
had  to  provide  a  clear  view  of  the  interior  of  the  combustor  over  all 
wavelengths  pertinent  to  thermal  radiation  measurements;  (2)  they  had  to 
withstand  the  thermomechanical  stresses  imposed  upon  them  by  the  high- 
pressure,  hot,  reactive  flame;  (3)  they  had  to  be  easily  removable  to  allow  for 
periodic  cleaning. 

Calcium  fluoride  was  chosen  as  the  window  material,  in  part  for  its  high, 
unvarying  transmissivity  between  0.4  and  8.0  /im.  The  transmissivity,  for  a 
10mm  thick  window,  is  approximately  85%  over  the  entire  wavelength  range 
[Lasers  and  App.,  1985].  The  actual  windows,  being  7.0  mm  thick,  have 
slightly  higher  transmissivities  than  the  10  mm  value.  However,  the  value  of 


Window 


Figure  5.  High-Pressure  Gaseous  Fuels  Combustor 
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transmissivity  taken  for  corrections  to  radiative  readings  was  that  of  the  10 
mm  windows. 
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The  second  optical  design  consideration  was  an  appropriate  window 
aperture— one  that  would  provide  observation  of  the  flame  through  a  known 
area  while  masking  the  hot  combustor  walls  from  the  detector  field  of  view. 
Obviously,  since  the  combustor  walls  would  be  hot  in  a  burning  environment, 
they  would  emit  radiation  in  the  same  spectral  region  as  the  flame  itself, 
interfering  with  attempts  to  determine  flame  radiation.  Figure  6  depicts  the 
geometry  discussed  and  gives  appropriate  dimensions. 


Window  overheating  was  avoided  through  the  use  of  water  cooling  and  a 
continuous  nitrogen  purge.  The  water  cooling  system  was  the  same  system' 
used  to  cool  the  combustor.  As  Figure  6  shows,  each  window  was  mounted  ■ 
on  a  stainless  steel  supporting  arm.  Combustor  cooling  water  flowed  through 
the  surrounding  water  jacket  cooling  the  window  support  arms  and  hence  the 
windows.  A  chromel/alumel  thermocouple  in  each  window  arm  provided  a 
continuous  readout  of  the  temperature  within  each  arm  to  ensure  that  it 
never  rose  above  35  C.  A  flow  of  dry  nitrogen,  injected  into  each  window 
arm  and  exhausted  to  the  combustor,  kept  soot,  water  vapor,  carbon  dioxide 
and  other  flame  gases  away  from  the  window  surfaces.  The  appropriate 
nitrogen  flow  rate  (0.05-0.10  g/s)  was  determined  by  adjusting  the  flow 
control  valve  until  the  radiation  reading  was  at  its  peak.  Observation  of  the 
flame  revealed  that  the  nitrogen  injection  only  slightly  altered  the  flow  field 
within  the  combustor. 

The  final  combustor  feature  is  the  system  controlling  combustion  pressure. 
It  consisted  of  a  back-pressure  valve  and  an  exhaust  bleed.  The  extent  to 
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which  the  valve  was  closed  and  the  mass  flow  rate  of  air  delivered 


determined  the  pressure  within  the  combustor.  Initial  high  pressure  runs 


revealed  difficulty  in  sustaining  a  flame  when  the  back-pressure  valve  came 


close  to  shutting  off  the  flow  entirely.  The  solution  was  to  include  an  exhaust 


bleed,  vented  directly  to  the  atmosphere.  Consequently,  a  1.27  cm  tube 


continuously  exhausted  a  portion  of  the  quenched  combustion  products 


thereby  reducing  the  sensitivity  of  combustor  performance  to  back-pressure 


valve  actuation. 


Igniter 


Ignition  was  accomplished  with  an  igniter  spark  arcing  across  a  3  mm  gap 
at  the  wall  of  the  combustor.  The  igniter  itself  consisted  of  a  tungsten- 


cathode  centered  in  a  stainless  steel  anode  with  the  two  electrodes  being 
insulated  from  each  other  by  ceramic  and  teflon  spacers.  A  110  V 


transformer  provided  a  5000  V  spark,  enough  to  jump  the  3  mm  gap  at 


pressures  as  high  as  0.5  MPa. 


Flame  ignition  was  assisted  by  heating  the  incoming  air  to  a  temperature 


of  500  K  or  higher.  The  ignition  procedure  consisted  of  warming  the  inlet  air 


to  at  least  500  K,  adding  fuel  at  stoichiometric  proportions,  positioning  the 


flameholder  immediately  upstream  of  the  spark  gap,  and  supplying  a  spark 


to  the  igniter.  A  stabilized  flame  was  detected  by  observing  a  rapid  increase 


in  the  output  of  a  radiometer  used  to  monitor  flame  zone  emission. 


Flameholder  and  Fuel  Injector 


If  the  flame  speed  is  less  than  the  flow  speed  of  fuel  and  air  in  a  subsonic 


flame  zone,  special  attention  must  be  given  to  keeping  the  flame  burning 
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stably.  In  these  experiments,  the  flow  velocity  was  set  at  6.0  m/s,  far  above 
the  laminar  flame  speed  of  0.4-0. 5  m/s.  This  combustor,  therefore,  employed 
a  flameholder  of  substantial  blow-off  speed-greater  than  25  m/s  at  standard 
temperature  and  pressure. 

The  flameholder  consisted  of  a  Hastalloy  disc  with  a  diameter  just  smaller 
than  the  inside  diameter  of  the  combustor  (3.8  cm).  The  disc  had  a  hoi"  of 
1.1  cm  diameter  at  its  center  through  which  the  fuel/air  mixture  entered  the 
combustion  zone.  See  Figure  7.  The  flameholder  could  be  translated  relative 
to  the  fuel  injector  to  provide  varying  degrees  of  premixing  between  the  fuel 


and  air. 
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Figure  7.  Fuel  Injector/Flameholder. 

The  fuel  injector  forced  the  fuel  and  air  to  mix  in  a  region  of  relatively 
high  velocity  immediately  upstream  of  the  flame  zone.  Since  the  speed  in  the 
mixing  region  was  approximately  fifteen  times  the  nominal  air  speed  in  the 


combustor,  the  danger  of  flashback  was  minimal. 


J. 


Because  the  flow  within  the  combustor  was  turbulent,  the  length  required 
for  fuel-air  mixing  was  expected  to  be  negligible.  Nevertheless,  experiments 
were  performed  to  verify  that  premixed  fuel  and  air  did  enter  the  flame  zone. 
To  accomplish  this,  a  stream  of  warm  nitrogen  was  fed  through  the  fuel 
injector  to  simulate  fuel  flow  while  ambient  air  flowed  through  the 
combustor.  The  difference  in  temperature  between  the  two  gases  was  20  K, 
easily  detectable  by  a  chromel/alumel  thermocouple.  Such  a  thermocouple 
was  translated  across  the  injector  outlet,  reading  temperatures  at  15  different 
radial  locations.  For  all  equivalence  ratios  from  0.5  to  2.0  and  for  nominal  air 
speeds  of  5,  10,  and  15  m/s  the  temperature  profiles  were  uniform  to  within 
1.5%,  indicating  that  the  two  gases  were  fully  premixed. 

Originally,  both  diffusion  and  premixed  flames  were  to  be  investigated: 
with  the  fuel  injector  designed  to  make  this  possible.  Early  runs  showed, 
however,  that  soot  was  present  at  nearly  all  conditions  investigated  in 
diffusion  flames,  so  this  portion  of  the  program  was  not  pursued. 

Flame  Positioning 

One  goal  of  this  research  was  to  study  soot  formation  and  flame  radiation 
as  a  function  of  distance  downstream  of  the  flameholder.  One  way  to 
provide  radiometer  access  to  several  flame  locations  would  be  to  build  several 
pairs  of  observation  ports.  Many  ports  in  the  combustor  wall  would  disturb 
the  flow,  however,  and  some  arrangement  for  moving  the  radiometer’s  field  of 
view  (e.g.  mirrors)  would  have  to  be  developed.  Moreover,  the  device  used  to 
monitor  flame  temperature  would  have  to  be  relocated  in  order  to  measure 
the  temperature  at  the  axial  location  where  radiometric  data  were  being 
taken.  The  solution  was  to  use  a  single  pair  of  windows,  a  fixed 


thermornetric  probe  location  and  a  sliding  flameholder/fuel  injector.  In  this 
way  the  flame  could  be  scanned  continuously  for  soot,  radiation,  and 
temperature,  all  from  a  single  vantage  point.  Flow  disturbances  resulting 
from  multiple  windows  were  also  minimized.  A  further  advantage  of  the 
sliding  flame  was  that  the  peak  of  radiation  would  never  be  lost  due  to  its 
occurrence  at  a  location  between  windows. 

Flame  positioning  hardware  consisted  of  a  gas-over-hydraulic  system, 
piping  and  solenoid  values.  Power  for  the  slide  arrangement  was  furnished 
by  a  hydraulic  cylinder,  with  the  hydraulic  fluid  pressurized  by  dry  nitrogen. 
The  liquid-gas  interfaces  lay  inside  two  pipes,  mounted  vertically. 
Pressurization  of  one  pipe  forced  hydraulic  fluid  into  one  end  of  the  cylinder,, 
driving  the  slide  in  one  direction.  Likewise,  pressurization  of  the  other  pipe 
drove  the  slide  in  the  opposite  direction. 

Flame  position  was  monitored  remotely  on  a  4-20  mA  meter  located  in  the 
control  room.  The  meter  was  connected  to  a  potentiometer,  mounted  next  to 
the  slide  and  attached  to  a  pulley  as  shown  in  Figure  8. 


Flame  Location  Indicator 


Figure  8.  Hydraulic  Slide  Mechanism. 


MEASUREMENTS  AND  DIAGNOSTICS 


Descriptions  of  all  devices  used  to  take  data  during  an  experimental  run 
are  presented  in  this  section.  Calibration  procedures  are  also  described  for 
each  device. 

Flame  Temperature 

Flame  temperature  is  important  when  analyzing  both  soot  thresholds  and 
flame  radiation.  Soot  thresholds  depend  exponentially  on  flame  temperature, 
through  the  rates  of  soot  precursor  formation  and  destruction.  Flame ' 
radiation,  in  contrast,  is  proportional  to  T4  for  a  gray  medium  and  may  also 
exhibit  a  weaker  dependence  since  gas  emissivities  depend  on  temperature. 

Two  methods  were  originally  considered  for  measuring  flame 
temperatures.  Both  were  line-of-sight  optical  techniques.  One,  the  Schmidt 
technique,  is  useful  for  isothermal  gray  gases.  The  second,  the  monochromatic 
emissivity  technique,  is  applicable  to  a  sufficiently  thick  chamber  of  gases 
containing  carbon  dioxide.  Brief  explanations  of  each  method  and  reasons 
for  rejecting  them  are  presented  below. 

If  a  gas/particle  mixture  is  gray,  i.e.  has  equal  emissivities  and 
absorptivities  at  all  wavelengths,  then  the  value  of  its  overall  emissivity  can 
be  determined  by  an  absorption  measurement.  Three  radiation  readings  are 
required: 
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(1)  emission  from  the  flame 

(2)  emission  from  a  com  ation  of  the  flame  and  a  blackbody  source 
placed  behind  the  flame 

(3)  emission  from  the  blackbody  source 


A  simple  calculation  then  yields  flame  emissivity.  Once  the  emissivity  is 
known,  it  can  be  used  to  calculate  the  flame  temperature  [Daws  and  Thring, 
1955].  This  procedure,  known  as  the  Schmidt  technique,  was  rejected 
because  it  requires  a  gray  medium.  This  may  be  an  acceptable 
approximation  for  very  sooty  flames,  but  in  this  work  most  of  the  flames  were 
either  completely  soot-free  or  at  the  verge  of  soot  formation;  therefore,  the 
gray  assumption  does  not  apply  for  most  flames  and  the  resulting 
temperature  calculations  would  have  been  erroneous. 


The  monochromatic  emissivity  technique  uses  an  optical  filter  or 
monochromator  to  isolate  radiant  energy  at  a  discrete  wavelength.  If  the 
flame  gas  is  known  to  have  an  emissivity  of  unity  at  that  wavelength 
measurement  of  the  monochromatic  emissive  power  allows  calculation  of  the 
flame  temperature  [Claus,  1981].  Carbon  dioxide  is  known  to  have  a  very 
strong  absorption  band  at  4.3  so  only  a  short  path  length  is  required  to 
ensure  an  emissivity  of  unity.  Such  conditions  occur  for  heavy  hydrocarbon 
fuels,  but  are  not  always  observed  when  burning  lighter,  gaseous  fuels.  The 
monochromatic  emissivity  technique  was  rejected  for  that  reason. 


Obviously,  neither  of  the  techniques  originally  proposed  was  applicable  to 
this  study.  This  was  due,  in  part,  to  a  shift  in  program  goals  from  use  of 
liquid  fuels  to  use  of  gaseous  fuels.  Consequently,  an  alternative  temperature 
measurement  technique  was  employed:  a  pyrometric  probe.  The  probe, 


termed  a  suction  thermometer,  samples  name  gases  through  a  choked-flow 
orifice  at  the  centerline  of  the  combustor  and  passes  them  through  a 
rotameter,  yielding  a  temperature  measurement  via  a  mass-flow 
measurement.  The  orifice  is  located  in  the  tip  of  a  water-cooled,  stainless 
steel  probe,  which  is  inserted  through  the  back  of  the  combustor  with  its 
inlet  immediately  downstream  of  the  observation  ports.  Temperature  is 
determined  from  a  mass  flow  measurement  by  using  the  relationship  between 
upstream  stagnation  pressure  (combustor  pressure  in  this  case)  and  mass  flow 
through  the  choked  orifice.  See  Figure  9. 
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Figure  9.  Suction  Thermometer  Inlet. 

The  probe’s  nozzle  was  carefully  shaped  in  an  effort  to  achieve  frictionless 
flow  in  the  converging  section  and  sonic  flow  within  the  throat  (diameter  = 


0.0838  cm).  Standard  construction  codes  [ASME,  1971]  call  for  a  circular 
profile  for  the  nozzle  inlet  when  sectioned  along  the  flow  axis  and  viewed 
from  the  side.  This  contour  was  approximated  by  a  rounded  inlet,  a  straight 
converging  taper,  and  a  rounded  transition  from  taper  to  constant-area 
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throat.  The  probe  also  incorporated  a  sudden  expansion  downstream  of  the 
throat  in  order  to  maintain  sonic  flow  in  the  constant-area  throat. 

Temperatures  could  be  measured  at  any  flame  location  visible  through  the 
window  because  the  probe  inlet  was  positioned  immediately  downstream  of 
the  observation  plane.  Axial  flame-temperature  profiles  were  collected  by 
traversing  the  flame  back  and  forth  relative  to  the  (fixed)  probe.  This 
accomplished  two  objectives.  First,  the  temperature  of  the  sooting  region 
could  be  examined  for  its  effects  on  soot  thresholds.  Second,  overall 
emissivities  could  be  mapped  out  from  a  knowledge  of  flame  temperature  and 
flame  radiation. 

A  test  was  performed  to  determine  how  probe  temperature  measurements 
compared  with  theory.  Air  of  known  temperature  and  pressure  was  sampled 
by  the  probe  with  the  volume  flow  rate  measured  using  the  rotameter.  The 
mass  flow  rate  was  compared  to  P0-T~1//2,  as  suggested  by  a  one-dimensional 
gas  dynamic  analysis.  Although  the  flow  rate  was  lower  than  predicted,  it 
was  proportional  to  P0*TJ"1//2  for  pressures  between  0.10  and  0.81  MPa  and 
temperatures  between  300  and  600  K. 

A  second  test  was  performed  to  avoid  extrapolating  low  temperature 
calibration  results  to  high  temperatures.  This  was  accomplished  by 
reconfiguring  the  rig  in  such  a  way  that  flame  temperatures  could  also  be 
monitored  by  a  Pt/Pt-10%Rh  thermocouple  inserted  into  the  combustion 
zone  through  an  observation  port.  A  lean  flame  was  established,  and  the 
thermocouple  reading  was  compared  to  simultaneous  suction  thermometer 
results.  Thermocouple  data  were  converted  to  temperatures  and  corrected 
for  radiative  losses.  Because  the  calibration  flames  were  lean,  and  therefore 


non-sooting,  the  emissivities  of  the  thermocouple  junctions  were  not  altered 
by  soot  deposition.  As  a  further  precaution,  the  thermocouple  was  coated 
with  a  thin  layer  of  Si02  to  prevent  catalytic  reactions.  A  calibration  curve 
resulted  from  these  tests,  allowing  determination  of  temperatures  under 
sooting  conditions  from  a  knowledge  of  the  probe  mass  flow  rate. 

Calculated  Flame  Properties 

Some  properties  used  in  the  analysis  of  the  flame  data  were  calculated 
rather  than  measured.  For  instance,  the  estimation  of  flame  emissivity  by 
Hottel’s  method  required  a  knowledge  of  partial  pressures  of  carbon  dioxide 
and  water  vapor,  while  analysis  of  the  critical  equivalence  ratio  data  required 
values  for  the  hydroxyl  radical  mole  fraction.  These  values  were  calculated 
for  equilibrium  conditions  at  the  measured  flame  temperature  and  combustor 
pressure.  In  addition,  the  adiabatic  flame  temperature  and  equilibrium 
species  mole  fractions  at  that  temperature  were  calculated  for  use  in 
comparisons  with  the  values  at  measured  flame  temperatures.  All  quantities 
were  determined  using  a  thermodynamic  equilibrium  computer  code  [Gordon 
and  McBride,  1973]. 

Air  Mass  Flow 

Combustion  air  flow  was  measured  using  an  8.86  mm,  0  =  0.2  orifice 
plate  in  conjunction  with  a  Vaiidyne  model  DP15-50  variable  inductance 
pressure  transducer.  See  Figure  10.  The  flow  through  the  orifice  plate  was 
assumed  to  be  incompressible  and  single-phase.  The  temperature  at  the 
orifice  location  was  experimentally  verified  to  be  the  same  from  day  to  day. 


Figure  10.  Fluid  Measurement  Orifice. 


The  associated  signal  analyzer,  a  Validyne  model  CD280,  delivered  a  DC 
voltage  in  response  to  a  given  pressure  drop  across  the  orifice  plate.  The 
following  equations,  taken  from  Goldstein  [1983],  give  the  flow  in  g/s  when  p 
is  in  psi  and  Ap  is  in  inches  of  water/15: 


m  =  pV 

V  =  a-cd 


2  Ap 
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1/2 


where 

m  =  mass  flow  rate 
v  =  volumetric  flow  rate 


A  =  area  at  orifice  (  =  —  d2) 

v  4  ' 

Cj)  =  discharge  coefficient  of  the  orifice 


(4) 

(5) 


•IvTOI 


amagaa 


Ap  =  pressure  drop  across  orifice 


p  =  gas  density 


/ 3  =  ratio  of  orifice  diameter  to  internal  pipe  diameter 


Equations  (4)  and  (5)  give 


m  =  c(p‘Ap)1//2 


where 


c  =  encompasses  all  constants  in  (5) 


p  =  pressure  upstream  of  orifice 


The  peculiar  units  for  Ap  are  a  matter  of  experimental  convenience. 


Calibration  with  a  mass  flow  sensor  (Micro  Motion  model  D,  accurate  to 
0.4%  of  reading)  showed  that  Eqn  (6)  gave  accurate  values  for  the  air  flow 


Combustor  Pressure 


Provision  was  made  to  remotely  monitor  the  combustor  pressure  by  using 


an  electromechanical  transducer.  A  pressure  tap  was  drilled  normal  to  the 


combustor  wall  with  the  pressure  sampled  through  a  6-mm  stainless  steel 


tube  capped  by  a  0-172  MPa  pressure  transducer  (Sensotec  model  A-10).  An 
electrical  signal  ran  from  the  transducer  to  a  digital  voltmeter  in  the  control 


room.  Calibration  was  accomplished  by  comparing  readings  of  combustor 


pressure  from  a  previously  calibrated  dial  gauge  with  the  voltage  output  of 


the  transducer. 


Combustor  Temperature 


Chromel/alumel  thermocouples  (Omega  Type  K)  were  used  for  all 
temperature  measurements,  except  flame  temperatures.  The  temperature 
was  monitored  at  two  locations  within  the  combustor:  in  an  insulated  section 
of  pipe  ten  diameters  upstream  of  the  water  jacket  and  in  an  uninsulated 
section  about  one  diameter  upstream  of  the  water  jacket.  The  former 
recorded  combustor  inlet  temperature.  The  latter  was  used  during  light-up 
to  indicate  when  a  flame  stabilized  on  the  flameholder  and,  because  of  its 
proximity  to  the  burning  zone,  could  also  reveal  flashback. 

Fuel  Mass  Flow 

Fuel  flow  rate  measurements  were  also  made  using  an  orifice  plate.  The' 
orifice,  with  a  pressure  tap  one  pipe  diameter  upstream  and  another  tap 
one-half  diameter  downstream,  was  installed  between  the  fuel  heater  and  the 
fuel  on/off  valve  (see  Figure  3).  The  fuel  orifice  dimensions  were  smaller  than 
those  of  the  air  orifice  to  allow  for  precise  measurement  of  the  small  flow 
rates  involved.  (The  inner  diameter  of  the  tube  was  1.09  cm  and  the  orifice 
diameter  was,  0.45  cm.)  The  pressure  upstream  of  the  orifice  and  the 
pressure  drop  across  it  were  monitored  by  Validyne  pressure  transducers  with 
the  resulting  signals  read  from  voltmeters. 

The  fuel  flow  rate  was  calibrated  as  a  function  of  upstream  pressure  and 
pressure  drop  across  the  orifice.  To  calibrate  the  pressure  upstream  of  the 
orifice  a  known  pressure  of  nitrogen  was  admitted  to  the  fuel  line  and  the 
flow  control  valve  was  shut  so  that  a  steady  stagnation  pressure  could  be 
read  by  the  upstream  transducer.  Its  output  voltage  was  then  read  from  a 


voltmeter. 


The  flow  rate  calibration  employed  both  the  upstream  pressure  and 
pressure  drop  transducers.  A  tank  of  known  volume  was  filled  with  nitrogen 
from  the  fuel  delivery  system  purge  line.  As  the  nitrogen  flowed  into  the 
tank,  the  pressure  and  pressure  drop  across  the  orifice  were  noted.  The 
gauge  pressure  after  filling  the  tank  and  the  time  to  fill  were  also  noted.  The 
mass  of  nitrogen  in  the  tank  was  then  calculated  as  follows: 


m  = 


Pg'V 

RT 


(7) 


where  m  =  mass  of  nitrogen  in  the  tank 
pg  =  gauge  pressure  after  a  trial 
V  —  tank  volume 
R  =  gas  constant  for  nitrogen 
T  =  absolute  temperature  of  the  nitrogen 
The  resulting  mass  was  divided  by  the  tank  fill  time,  yielding  the  mass  flow 
rate.  The  mass  flow  rates  of  all  trials  were  plotted  against  the  square  root  of 
the  product  of  upstream  pressure  and  orifice  pressure  drop: 


nitrogen  =  019  (P’Ap)1/2  (8) 

where 

^nitrogen  ~  ma8S  flow  rate  of  nitrogen  (g/s) 
p  =  upstream  pressure  readout  (V) 

Ap  =  pressure  drop  readout  (mV) 

Solving  for  Ap  gave 


Ap 


27.7  m 


2 

nitrogen 


(9) 


P 


43 


r  J 

% 


The  desired  value  of  m  was  determined  from  the  nominal  combustor  air 
flow  rate  and  desired  equivalence  ratio.  The  readout  of  upstream  pressure 
was  then  used  to  give  the  necessary  pressure  drop,  which  was  scaled  to 
account  for  differences  between  the  molecular  weight  of  nitrogen  and  that  of 
the  fuel  to  be  supplied.  The  procedure  is  illustrated  by  the  following 
expressions. 


nVuel  mnitrogen 


^fuel 


(  ^nitrogen 


1/2 


(10) 


ramtrogen 


M, 


T 

fuel  nitrogen 


M 


nitrogen 


L  fuel 


1/2 


(11) 


for  the  same  upstream  pressure  where 
m  =  mass  flow  rate 
p  —  gas  density 

M  =  molecular  weight  of  the  gas 
T  =  gas  temperature 

The  nitrogen  was  at  a  temperature  of  300  K  while  the  fuel  was  held  at  375  K 
by  means  of  the  fuel  heater.  Scaling  of  the  temperatures  then  yielded: 


ihfue,  =  0-17  (p-Ap)1/2  (12) 

Radiometer 

Proper  measurement  of  flame  radiation  called  for  careful  selection  of 
radiometric  components.  The  four  principal  components  were  the  power 
meter,  the  rotating  wheel  chopper,  the  combustor  window  assemblies,  and  the 
calibration  furnace.  The  power  meter  was  a  Laser  Precision  model  RL-3610 
of  1.0  cm2  collection  area  used  in  conjunction  with  a  Laser  Precision  model 
CTX-530  chopper  operating  at  30  Hz.  It  was  able  to  detect  radiant  energy 


from  0*4  ^ m  to  12.0  fim.  This  wide  response  range  is  useful  only  if  the 
window  separating  the  power  meter  and  the  interior  of  the  combustor  is 
transparent  to  radiation  throughout  the  same  wavelength  band. 
Consequently,  calcium  fluoride  was  chosen  for  the  window  material  since  it 
transmits  approximately  95%  of  the  incident  energy  from  0.7  /im  to  8.0  /im. 
One  important  consideration  concerns  the  space  between  the  window  surface 
and  the  flame  itself  (Figure  6).  Had  this  space  been  allowed  to  fill  with 
carbon  dioxide  or  water  vapor  from  either  the  atmosphere  or  the  flame, 
radiative  emission  from  the  flame  would  have  been  attenuated.  Because  this 
path  was  continuously  purged  with  nitrogen,  however,  no  such  attenuation 
occurred  (nitrogen  being  transparent  to  infrared  radiation).  The  final 
component  of  the  radiometric  system  was  an  Infrared  Industries  model  IR- 
464  blackbody  furnace  with  a  model  IR-101C  controller.  With  a  maximum 
cavity  temperature  of  1473  K  this  blackbody  provided  a  standard  for 
calibration  of  the  radiation  power  meter. 

Calibration  was  accomplished  by  placing  the  radiometer  and  chopper 
directly  in  front  of  the  blackbody  furnace  at  a  known  distance  of  separation 
and  measuring  the  output  voltage  of  the  radiometer  at  various  blackbody 
temperatures.  The  following  analysis  describes  the  process  in  detail. 


<k-r  =  W^e-r 

=  I.  A.  A,  /  il. 


where 


A  =  area  of  emitting  or  receiving  surface 
d  =  distance  from  emitter  to  receiver 


I  =  normal  radiative  intensity  from  the  emitting  surface 
^e-r  =  solid  angle  between  emitter  and  receiver 
a  =  5.67  -1(T8  W/  m2-K4 
q  =  radiant  energy  emission  rate 
S  =  signal  from  the  radiometer 

and  the  subscripts  are 
e  =  emitter 
r  =  receiver 


For  a  blackbody 


where 

E  =  emitted  flux 

with  subscript 

b  =  blackbody 

Using  Eqns  (15)  and  (16) 

„  fTb*  Ai-A, 

bb  ”  '  > 2 

n  db-r 

For  the  flame  as  the  emitter  we  can  rearrange  Eqn  (14)  to  yield 

_  qf-r’df-r 
A,A, 

and  use  Eqns  (13)  and  (15)  to  give 

Sf 

*lf-r  c  flb-r 
bb 


RESULTS  AND  ANALYSIS 


This  portion  of  the  report  consists  of  two  sections.  The  first  is  devoted  to 
soot  threshold  equivalence  ratio  data  and  analysis.  The  second  discusses 
flame  radiation  results  and  analysis. 

Soot  Thresholds 

Figures  11  through  17  show  the  results  of  variations  in  pressure  on  the 
critical  equivalence  ratio,  4>c,  for  all  fuels  studied.  Critical  equivalence  ratio 
data  are  presented  at  multiple  temperatures  for  three  of  the  fuels,  ethane, 
propane,  and  ethylene.  Clearly,  both  combustor  pressure  and  inlet 
temperature  strongly  affect  0C. 

Effect  of  Pressure 

Examination  of  the  data  for  any  fuel  shows  that  the  critical  equivalence 
ratio  falls  as  combustor  pressure  rises.  This  agrees  with  the  work  of 
Fenimore  et  al.  [1957],  who  studied  the  subatmospheric  pressure  dependence 
of  4>c.  Further  examination  of  the  figures  reveals  that  at  higher  pressures,  0.4 
MPa  and  above  for  most  fuels,  the  critical  equivalence  ratio  becomes  less 
sensitive  to  pressure,  approaching  a  horizontal  asymptote.  For  methane, 
ethane,  ethylene,  and  propylene  the  high-pressure  limiting  value  of  0C  is 
approximately  1.20,  when  the  inlet  air  temperature  is  700  K. 
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Figure  11.  Methane:  Critical  Equivalence  Ratio  vs.  Pressure. 
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Figure  12.  Ethane:  Critical  Equivalence  Ratio  vs.  Pressure 
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Figure  14.  Ethylene:  Critical  Equivalence  Ratio  V9.  Pressure 
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Figure  16.  C-4  Fuels:  Critical  Equivalence  Ratio  vs.  Pressure. 


For  propane  the  value  is  about  1.15,  for  cyclopropane  1.05,  and  for  the  four- 
carbon  fuels  it  appears  to  be  lower,  but  no  high-pressure  data  were  taken  to 
confirm  this.  The  fact  that  methane,  ethane,  and  propylene  all  approach  the 
same  value  of  d>c  at  high  pressures  is  interesting  because  these  fuels  differ 
markedly  in  their  sooting  tendencies  at  low  pressures.  Notice  that  methane 
did  not  produce  soot  at  any  equivalence  ratio  below  a  pressure  of  0.4  MPa. 
Similarly  ethane  produced  no  soot  below  a  pressure  of  0.3  MPa.  The  two 
olefins,  ethylene  and  propylene,  both  formed  soot  over  the  entire  pressure 
range  (0.1  to  0.8  MPa). 

Effect  of  Temperature 

Three  fuels,  ethane,  propane,  and  ethylene,  were  burned  at  different  inlet 
air  temperatures,  with  all  other  fuels  run  at  an  inlet  temperature  of  700  K 
only.  As  Figures  12,  13,  and  14  show,  a  change  in  inlet  air  temperature  can 
change  the  value  of  4>c  at  a  given  pressure.  The  <f>c  versus  pressure  behavior 
is,  however,  qualitatively  the  same  for  all  inlet  temperatures.  An  interesting 
feature  of  the  temperature  curves  is  that  for  different  fuels  the  inlet 
temperature  at  which  0C  is  a  minimum  varies.  Also,  no  simple  relationship 
exists  relating  <f>c  to  inlet  air  temperature  alone.  Propane,  for  example,  shows 
a  decrease  in  soot  threshold  as  the  inlet  air  temperature  rises  to  600  K.  The 
curve  at  700  K,  however,  lies  above  the  lower  temperature  curves.  Ethane, 
on  the  other  hand,  has  its  lowest  4>c  values  at  an  inlet  temperature  of  500  K, 
while  ethylene’s  lowest  <t>c  values  are  for  inlet  temperatures  of  500  and  600  K. 
Our  results  for  ethylene,  the  only  fuel  studied  as  such  in  another  work,  show 
the  same  behavior  as  reported  by  Wright  [1969]. 


A /  A J1,  A/,  V,  /.  WkA  if  .V-/* 


Effect  of  Fuel  Type 


As  mentioned  above,  the  various  fuels  had  widely  varying  soot  thresholds. 
For  the  purposes  of  comparison  we  will  look  at  the  different  results  for  inlet 
temperatures  of  700  K.  More  general  comparisons  will  be  made  later  in  the 
analysis  section  of  this  chapter. 

Several  fuel  parameters  have  been  suggested  as  affecting  critical 
equivalence  ratio.  For  instance,  a  comparison  of  the  <f>c  versus  pressure 
curves  for  the  parafinic  fuels  (methane,  ethane,  propane,  and  n-butane) 
shows  an  increasing  tendency  to  form  soot  as  molecular  weight  increases 
(from  16  for  methane  to  58  for  n-butane).  Similar  behavior  is  observed  for  the 
olefins  (ethylene,  propylene,  and  1-butylene)  from  ethylene  up  to  1-butylene'.' 
A  comparison  of  fuels  with  the  same  number  of  carbon  atoms  shows  that 
paraffins  generally  soot  less  readily  at  low  pressures  than  olefins,  a  result 
contrary  to  many  atmospheric  pressure  studies  (reviewed  in  Calcote  and 
Manos  [1983]).  For  the  three-carbon  fuels  we  can  make  the  further 
comparison  with  cyclopropane,  showing  that  cyclopropane  forms  soot  more 
readily  than  either  propylene  or  propane.  No  simple  relationship  is  apparent, 
however,  between  4>c  and  either  the  number  of  carbon-carbon  bonds  or  the 
H/C  ratio.  In  the  following  section  we  will  see  that  explicit  fuel  effects  are 
negligible  and  that  it  is  the  temperature  and  pressure  that  are  needed  to 
predict  4>c.  It  should  be  noted,  however,  that  fuel  type  does  have  an  implicit 
effect  on  <t>c  through  its  determination  of  the  equilibrium  hydroxyl  radical 
concentration. 


Soot  Threshold  Analysis 


Since  soot  threshold  data  depend  on  pressure  and  temperature,  any  model 
explaining  the  data  must  incorporate  these  variables.  Also,  in  accordance 
with  the  work  of  Harris  et  al.  [1986],  possible  fuel  chemistry  effects  should  be 
accounted  for.  The  analysis  presented  below  follows  closely  the  work  of 
Harris  et  al.  and  Takahashi  and  Glassman  [1984].  It  is  based  on  three 
assumptions.  First,  that  soot  formation  in  premixed  flames  is  directly  related 
to  the  formation  of  precursors,  and  that  only  two  chemical  reactions  are 
necessary  to  describe  the  evolution  of  the  precursor  concentration:  (i) 
formation  of  precursors  via  fuel  pyrolysis;  (ii)  oxidation  of  precursors  via 
hydroxyl  radicals.  Second,  that  the  hydroxyl  radical  concentration  can  be. 
described  by  its  equilibrium  value,  calculated  using  the  temperature 
measured  at  that  location  in  the  flame  where  soot  is  first  observed.  Third, 
that  the  concentration  of  fuel  is  linearly  proportional  to  the  equivalence 
ratio. 

Note  that  the  second  assumption  differs  from  the  works  of  Harris  et  al. 
and  Takahashi  and  Glassman.  Harris  et  al.  assume  <t>c  depends  on  the 
maximum  measured  flame  temperature  while  Takahashi  and  Glassman 
assume  xJj, c  depends  on  the  adiabatic  flame  temperature.  This  point  will  be 
discussed  again  in  a  subsequent  section. 

The  first  assumption  may  be  expressed  as 

[PI  =  kF  [F]  -  kD  [OH]  [P]  (21) 

where  [P]  is  the  time  rate  of  change  of  the  precursor  concentration,  [F],  [OH], 
and  [P]  are  the  fuel,  hydroxyl  radical  and  precursor  concentrations,  and  kF 
and  kD  the  rate  coefficients  for  precursor  formation  and  destruction.  Back 


reactions  have  been  neglected  in  both  cases.  Note  that  Eqn  (21)  is  based  on 
the  analysis  proposed  by  Harris  et  al.  and  Takahashi  and  Glassman,  which 
models  precursor  formation  as  a  first  order  reaction  and  precursor 
destruction  as  a  second  order  reaction. 

By  definition,  the  critical  equivalence  ratio  occurs  when  the  precursor 
concentration  is  at  steady  state.  Consequently 


[FU,  -  (OHl«,  [P]*,  7-  (22) 

kf 

The  second  assumption  allows  determination  of  [OH]0^,  using  a  standard 
chemical  equilibrium  code,  if  the  "correct"  temperature  is  known.  Several 
candidate  temperatures  have  been  suggested,  including  the  adiabatic  flame 
temperature  for  products  of  either  C02  and  HzO  [Olson  and  Pickens,  1984]  or 
CO  and  H20  [Millikan,  1962].  Another,  more  physical,  suggestion  is  an 
experimentally  measured  temperature  [Harris  et  al.;  MacFarlane  et  al.,  1964; 
Flower  and  Bowman,  1986]. 

The  final  step  in  the  analysis  is  to  relate  the  fuel  concentration  at  the 
critical  equivalence  ratio  to  the  critical  equivalence  ratio.  An  expression, 
similar  to  those  presented  by  Millikan,  Takahashi  and  Glassman  and  Harris 
et  al.,  is  derived  below. 

The  equivalence  ratio  is  commonly  defined  as 
mr/mi  ]  actual 

<t>  =  —r  , - -  (23) 

mp/mA]stoi 

where  mF  is  the  mass  of  fuel  present  and  mA  is  the  mass  of  air  present. 
Then, 
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^  _  XFMFMTmT/XAMAMTmT  !  act 
XFMA^TmT / XA^A^TmT  !  st°i 

with  Xf  an<^  Xa  tbe  m°le  fractions  of  fuel  and  air,  Mp,  MA  and  MT  the 
molecular  weights  of  the  fuel,  air  and  fuel/air  mix,  and  m?  being  the  total 
mass  (fuel  +  air).  Obviously,  Eqn  (24)  reduces  to 

4  =  (25) 

XF/XAIst01 

and  can  be  rewritten  as 


1-Xf 
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The  denominator  is  a  constant,  for  a  given  class  of  fuels.  For  an  alkane 
fuel,  for  instance,  having  C02  and  H20  as  products 

C„H,n+,  +  (^-1(0,  +  3.76  N,)  - 

nCOz  +  (n+l)H20  +  1.88(3n+l)N2  (27) 


Xf/Xa'i3^*  = 


l/(— (3n+l)  +  1) 

™  (3n+l)/(~~“(3n+l)  +  1) 


Equation  (28)  reduces  to 

Xf  ,  1 

-  stoi  =  - - - - 

XA  1  (3n+l)2.38 

If  il>  (the  equivalence  ratio  based  on  CO  and  H20)  is  used  instead, 


!  stoi  = 


(2n+l)2.38 
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,1 
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For  any  alkane,  Eqn  (26)  shows  Xf  <  0-2  if  4>  <  2.0  since  n  >  1.  The 
same  holds  true  for  1 1>  <  1.43.  For  acetylene  or  an  aromatic  fuel 


2n+  — 

CnHa  +  ( - — )(02  +  3.76  N2)  =  nC02  +  —  H20  +  4.7n  N2 

2  2 


] stoi  = 


^-(4.76) 

For  acetylene,  Eqn  (26)  shows  Xf  <  0.17  if  <f>  <  2.  Aromatic  fuels  would  have 
Xf  <  0.12  for  the  same  conditions.  Equation  (26)  then,  shows  that  in  general 


1-Xf 


where  c  is  a  constant  dependent  on  the  fuel.  Furthermore,  since  Xf  <  0.2, 
the  right  hand  side  of  Eqn  (33)  is  approximately  equal  to  Xf/c*  Thus, 

Xf  =  c  ^  (34) 

Notice  that  the  expressions  due  to  Takahashi  and  Glassman  and  Harris  et  al. 
implicitly  assume  constant  pressure  and  temperature  since  they  equate  the 
fuel  concentration,  instead  of  the  mole  fraction,  to  4>.  Their  assumption  of 
constant  temperature  and  pressure  is  reasonable  considering  their 
temperatures  varied  by  less  than  20%  and  all  their  results  were  for  a  single 
pressure.  In  contrast,  the  expression  given  in  Eqn  (34)  requires  neither 
constant  pressure  nor  temperature. 

Transforming  the  concentrations  in  Eqn  (22)  into  mole  fractions,  and 
recalling  that 


raw 


where  is  the  partial  pressure  of  species  i,  Ru  is  the  universal  gas  constant, 
T  is  the  temperature,  and  p  the  mixture  pressure  shows 


P  Xf  _  .  .  _  P  Xqh  P  Xp  kD 
RUT  “  1  ]  “  RUT  RUT  kF 

Using  Eqn  (34)  for  Xf>  an<^  dividing  both  sides  by  Xoh  an<*  c  yields 

4>c  _  kD  p  Xp 

Xoh  RuT  c 

The  rate  coefficients  are  then  written  in  Arrhenius  form 


(36) 


(37) 


kj  =  A;  e  EyRuT  (3*)-. 

where  A,  is  the  collision  factor  and  Ej  the  activation  energy.  Taking  the 
natural  logarithm  of  both  sides  gives 


_  Ef-Ed  Xp.  p  ad 
Xqh  RuT  n  c  RuT  Af 


(39) 


If  we  make  the  simplifying  assumption  that 


XpAp 

cAf 


is  independent  of 


temperature  and  pressure,  then  for  a  single  fuel 


,  ,  p  EF— Ep  Xp  Ap 

In  -  —  In  -f-  =  — — — - 1-  In - — — 

Xoh  T  RuT  c  AfRu 

Equation  (40)  forms  the  basis  for  analysis  of  our  critical  equivalence  ratio 

data.  In  what  follows  we  examine  the  outcome  of  the  application  of  Eqn  (40) 

to  the  critical  equivalence  ratio  data  presented  in  Figures  11  to  17,  using 

both  adiabatic  and  measured  flame  temperatures.  We  will  also  examine  the 

applicability  of  using  the  effective  equivalence  ratio,  0C,  as  proposed  by 

Takahashi  and  Glassman.  Note  that  neither  the  H/C  ratio  of  the  fuel,  nor 


any  other  fuel  property  is  included  in  the  <pc  figures.  This  is  in  contrast  to 
the  results  of  several  other  studies.  At  the  conclusion  of  this  section  we  will 
see  that  such  terms  are  insignificant  and  can  therefore  be  ignored. 

The  goals  in  applying  Eqn  (40)  are  to  show  that  the  two-step  model  for 
soot  precursor  evolution  correctly  predicts  soot  threshold  behavior  over  the 
pressure  range  0. 1-0.8  MPa,  and  to  determine  which  temperature  and 
equivalence  ratio  best  describe  the  data. 

Figure  18  shows  results  obtained  using  the  adiabatic  flame  temperature 
and  corresponding  equilibrium  OH  concentration  in  Eqn  (40).  The  curves  for 
the  various  fuels  are  not  linear  and  do  not  lie  on  top  of  one  another.  The  use 
of  the  measured  flame  temperatures  and  associated  equilibrium  OH  mo!e'; 
fractions,  however,  yields  very  good  agreement  between  Eqn  (40)  andi 
experimental  data.  This  is  illustrated  in  Figure  19.  A  linear  least  squares  fit 
of  these  results  gives 


=  0.040 
Ax 

where 


4  Tf 

y  =  In - 

Xoh  P 


The  correlation  coefficient,  r2,  is  0.967  for  measured  flame  temperature  data 
but  only  0.850  for  adiabatic  flame  temperature  data.  Thus,  the  measured 
flame  temperature  should  be  used  when  predicting  critical  equivalence  ratio 
behavior.  Significantly,  Harris  et  al.  arrived  at  a  similar  conclusion  regarding 
measured  versus  adiabatic  flame  temperatures,  although  they  compared  the 


y  - 


adiabatic  flame  temperature  with  the  maximum  measured  flame  temperature, 
instead  of  the  temperature  measured  at  the  point  of  incipient  sooting. 

It  can  easily  be  shown  that  the  results  presented  here  are  sufficiently 
sensitive  to  discriminate  between  the  use  of  measured  and  adiabatic  flame 
temperatures.  Starting  with  Eqn  (40)  and  transforming  it  to 


-  C  e-^/^  (41) 

XOH 

(where  C  is  a  temperature  independent  term  containing  the  pre-exponential 
factors  for  the  rates  of  precursor  formation  and  destruction,  the  precursor 
mole  fraction  and  the  constant  relating  equivalence  ratio  to  fuel 
concentration)  we  can  show  that  if  our  results  are  insensitive  to  the  choice  oT 
flame  temperature  either  value  used  in  the  argument  of  the  exponent  (Eqn 
(41))  would  yield  the  same  functional  form.  To  show  that  this  is  not  so  we 
substitute  the  measured  and  adiabatic  temperatures  into  Eqn  (41)  and 
equate  the  two  expressions.  After  canceling  all  temperature  independent 
terms,  we  see  that 

AE/RTmM,  __  AE/RT^ii  (42) 

Rewriting  Tadia 

Tadia  =  Tmeas  + 

and  substituting  into  Eqn  (42)  gives 
e-AE/RT^l-AT/Tldii) 

for  the  right  hand  side.  If  we  expand  the  quantity  (1  -  AT/Tadja)-1  and 
substitute  it  into  Eqn  (43),  we  get 


-AE/RT.*.  e— AEAT/RT ** 


w.v.v.v 
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If  the  measured  and  adiabatic  temperatures  give  equivalent  results,  we  can 
interchange  them.  The  left  hand  side  of  Eqn  (42)  will  cancel  with  the  first 
term  in  Eqn  (44),  leaving 

e— aeat/RT^.  =  l  (45) 

Consequently,  if  AT  is  quadratic  in  Tadia,  we  will  be  incapable  of 
distinguishing  between  the  two  temperatures.  Figure  20,  which  presents  AT 
versus  Tadia,  immediately  shows  that  AT  is  not  quadratic  in  Tadia  so  we 
conclude  that  our  results  are  sufficiently  sensitive  to  distinguish  between  the 
TmeM  and  Tadia. 

The  effective  equivalence  ratio  suggested  by  Takahashi  and  Glassman,  xJjc, 
was  also  employed  to  analyze  the  V>c  results.  Results  using  this  equivalence 
ratio,  considered  by  those  authors  to  be  more  pertinent  to  rich  flames,  are 
presented  in  Figures  21  and  22.  Figure  21  presents  results  using  measured 
flame  temperatures  while  Figure  22  presents  results  using  adiabatic  flame 
temperature.  The  correlation  coefficients  from  least  squares  linear  fits  are 
0.967  and  0.863,  respectively. 

From  this  analysis  we  can  see  that  the  use  of  Eqn  (40)  with  measured. 
flame  temperatures  and  either  V>c  or  V’c  yields  the  best  results.  Significantly, 
the  combination  of  either  4>c  or  Vfc  and  adiabatic  flame  temperature  is  shown 
to  give  poor  results. 

Several  groups  have  suggested  V>c  (or  VO  should  depend  explicitly  on  some 
"fuel  parameter,"  such  as  the  H/C  ratio  [Harris  et  al.]  or  the  number  of 
carbon-carbon  bonds  [Takahashi  and  Glassman].  Consequently,  Eqn  (40)  was 
modified  to  include  an  explicit  dependence  of  V>c  on  parameters  such  as  H/C 
ratio,  H/C  mass  ratio,  number  of  carbon-carbon  bonds,  number  of  carbon 
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Figure  20.  Temperature  Difference  vs.  Flame  Temperature,  Squared 


atoms  and  fuel  molecular  weight. 

The  data  were  reanalyzed  with  the  results  presented  in  Table  1  and 
Figure  23.  It  is  readily  apparent  though,  that  none  of  the  parameters 
exhibited  improved  agreement  between  the  experimental  results  and  Eqn  (40) 
for  the  fuels  and  conditions  studied  here.  We  thus  conclude  that  </>c  exhibits 
no  explicit  dependence  on  any  "fuel  parameter."  However,  fuel  type  does 
have  an  implicit  effect  on  d>c  through  its  impact  on  the  equilibrium  hydroxyl 
radical  concentration. 

In  summary  several  conclusions  are  apparent.  First,  the  measured  flame 
temperature  predicts  critical  equivalence  ratios  better  than  the  adiabatic 
flame  temperature.  In  addition,  the  measured  temperature  should  be  that 
temperature  occurring  at  the  location  where  soot  is  first  observed  to  form. 
This  result  i?  intuitively  obvious  since  it  should  be  the  temperature  occurring 
at  the  point  where  the  precursor  concentration  is  in  steady  state  that 
controls  the  chemical  kinetic  processes  leading  to  incipient  sooting.  Second, 
either  the  standard  equivalence  ratio,  4>,  or  the  equivalence  ratio  based  on 
CO  and  H20  as  products,  ipc,  accurately  predicts  incipient  sooting  behavior. 
Finally,  for  the  fuels  and  conditions  investigated  here  no  explicit  dependence 
of  <f)c  (or  rjjc)  on  any  "fuel  parameter"  is  apparent,  if  the  temperature 
occurring  at  the  point  of  incipient  soot  formation  is  used  in  Eqn  (40).  No  fuel 
parameter  dependence  inserted  into  a  modified  version  of  Eqn  (40)  gave  any 
improvement. 

The  conclusions  drawn  here  imply  that  the  mechanism  leading  to  soot 
formation  is  similar  for  all  fuels.  This  hypotheses  should  be  investigated 
further  to  confirm  its  validity. 
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Flame  Radiation 

Figures  24  through  27  depict  radiation  intensity  versus  axial  location  for 
ethane,  propane,  ethylene,  and  propylene  flames,  respectively.  The  inlet 
temperature  was  700  K  for  all  tests,  while  the  pressure  was  varied  from  0.4 
to  0.8  MPa  (except  for  propylene,  for  which  no  0.8  MPa  data  were  taken).  In 
addition,  tests  were  run  at  three  equivalence  ratios,  namely:  0.80,  1.00,  and 
1.25.  These  conditions  produced  flames  of  negligible  or  zero  soot.  Flame 
radiation  intensity  was  then  plotted  vs  x/D. 

To  obtain  mappings  of  radiation  intensity  versus  axial  location, 
measurements  were  taken  at  integral  values  of  x/D  in  most  cases;  however, 
at  each  pressure  the  radiometric  output  was  monitored  as  a  continuous 
function  of  axial  location  for  25%  of  the  cases  so  that  curves  could  be  fitted 
accurately  between  integral  x/D  locations  for  all  conditions. 

Three  quantities  were  varied  to  examine  their  effects  on  flame  radiation 
levels.  Two  of  these  were  combustor  pressure  and  equivalence  ratio,  while 
the  third,  fuel  type,  covered  two  paraffins,  ethane  and  propane,  and  two 
olefins,  ethylene  and  propylene.  Note  that  propane  and  ethylene  have  the 
same  number  of  carbon-to-carbon  bonds.  The  effects  of  the  above 
parameters  on  flame  radiation  intensity  are  discussed  below. 
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To  simplify  the  discussion  and  to  focus  on  the  most  important  data,  the 
analysis  is  restricted  to  the  peak  levels  of  radiation  for  each  axial  profile. 


Effect  of  Pressure 


An  increase  in  combustion  pressure  yielded  a  linear  increase  in  peak 
intensity  for  all  fuels,  as  shown  in  Figures  24  through  27.  However,  Figure  28 
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shows  that  the  slope  of  the  line  denoting  peak  intensity  versus  pressure 
differs  between  equivalence  ratios.  The  slope  is  1.0  for  an  equivalence  ratio  of 
unity,  falling  to  0.90  for  <t>  =  1.25  and  0.70  for  <f>  =  0.80.  Evidently,  the 
changes  in  flame  gas  density,  flame  temperature,  species  concentrations,  and 
emission  bandwidths  brought  about  by  changes  in  pressure  are  dependent  on 
the  equivalence  ratio  of  the  unburnt  mixture.  It  is  interesting  to  note  that 
despite  the  complicated  nature  of  the  pressure  effect,  the  intensity  responds 
linearly  to  pressure  changes  in  the  range  0.4  to  0.8  MPa.  This  linear 
response  agrees  with  results  taken  from  liquid  fuel  studies  for  the  pressure 
range  covered  here  [Claus;  Naegeli  et  al.,  1983].  It  is  also  corroborated  by 
the  results  of  MacFarlane,  who  burned  gaseous  methane  in  air. 


Effect  of  Equivalence  Ratio 


The  radiation  intensity  follows  several  clear  trends  as  equivalence  ratio 
varies.  First,  equivalence  ratios  of  0.80  and  1.25  have  nearly  equal  peak 
intensities.  This  is  shown  in  Figure  29,  which  depicts  the  maximum  radiant 
intensity  averaged  over  all  fuels  as  functions  of  4>  and  p.  In  general,  about 
20%  more  C02  and  20%  less  H20  is  produced  for  <2>  =  0.80  than  for  <j>  = 
1.25  at  equilibrium  flame  conditions,  as  determined  using  the  equilibrium 
code  of  Gordon  and  McBride.  The  temperatures  in  the  experiment  discussed 
here  were,  on  the  average,  10%  higher  for  6  =  1.25  than  for  <t>  =  0.80,  as 
demonstrated  in  Figures  30  and  31.  An  important  conclusion  is  that  the 
higher  C02  concentration  give  the  lean  flame  as  great  an  emissive  flux  as  the 
rich  flame,  despite  its  higher  H20  concentration  and  temperature.  This  is 
logical  because  the  infrared  emissivity  of  C02  is  much  higher  than  that  of 
H20  [Siegel  and  Howeilj.  Also,  as  shown  in  Figure  32,  the  stoichiometric 


Figure  24.  Ethane:  Flame  Radiation  Profiles. 
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Figure  25.  Propane:  Flame  Radiation  Profiles. 
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Figure  28.  Peak  Flame  Radiation  Intensity  vs.  Pressure. 
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flames  exhibited,  on  average,  a  20%  higher  intensity  than  flames  burning  at 
equivalence  ratios  of  0.8  and  1.25.  This  is  also  reasonable  because  the 
equilibrium  H20  concentrations  of  the  stoichiometric  flame  were  nearly  equal 
to  those  for  <f  =  1.25,  and  the  equilibrium  C02  concentrations  were  10% 
higher  than  the  values  for  &  =  0.80.  Moreover,  stoichiometric  flame 
temperatures  were  generally  5-10%  higher  than  for  4>  =  1.25.  One 

observation,  then,  is  that  equivalence  ratios  of  0.80  and  1.25  yield  nearly 
equal  radiation  intensities,  which  are  about  20%  lower  than  for 
stoichiometric  mixtures. 

Another  observation  is  that  flame  length  increases  with  equivalence  ratio, 
so  that  the  peak  intensity  moves  downstream  as  the  equivalence  ratio  is 
increased.  This  effect  can  be  seen  very  clearly  in  a  comparison  of  any  curves 
depicting  4>  =  0-80  and  <j>  =  1.25  (e.g.  Figure  25).  The  former  peak  early, 
£  usually  before  x/D  =  1.5;  the  f>  =  1.25  curves,  on  the  other  hand,  continue 

to  rise,  peaking  at  x/D  =  4.  The  explanation  is  that  rich  mixtures  require  a 
longer  time  {and  therefore  length)  to  achieve  complete  combustion  because 
more  fuel  is  burnt  per  unit  volume  of  air. 

Effect  of  Fuel  Properties 

The  highest  flame  radiation  levels  were  obtained  with  propylene,  followed 
by  propane,  ethylene,  and  then  ethane.  These  results  indicate  that  flame 
radiation  may  depend  on  the  hydrogen-to-carbon  (H/C)  ratio  of  the  fuel,  as 
suggested  by  gas  turbine  experiments  [Naegeli  et  al.,  1982;  1983],  or  on  the 
number  of  bonds  in  the  fuel  molecule.  Since  the  peak  radiation  levels 
appeared  to  be  a  function  of  the  hydrogen-to-carbon  ratio  of  the  fuel,  Ima, 
was  plotted  against  H/C.  As  Figure  33  shows,  the  radiant  intensity  generally 


falls  with  increasing  H/C,  which  agrees  with  the  results  of  Moses  and  Naegeli 
[1979).  This  agreement  indicates  that  although  gas  turbine  flame  radiation  is 
certainly  complicated  by  liquid  fuel  spray  phenomena,  it  is  also  dependent 
upon  the  degree  of  fuel  molecule  saturation. 

Because  the  curves  in  Figure  33  are  qualitatively  similar  but  vertically 
separated,  other  expressions  were  tested  to  determine  if  they  more  accurately 
described  the  data.  Recalling  that  density  played  a  role  in  flame  gas 
concentrations  for  the  soot  threshold  analysis,  it  might  be  expected  to 
influence  the  C02  and  H20  concentrations  in  the  same  manner.  A  linear 
dependence  of  radiation  intensity  on  gas  density  is  not  anticipated  because  of 
self  absorption  along  the  radiation  path  length,  but  a  gas  of  higher  density 
should  emit  more  energy  than  one  of  low  density.  Figure  34  shows  that  the 
families  of  curves  lie  slightly  closer  together  when  corrected  for  differences  in 
density,  indicating  that  the  gas  density  may  be  an  important  factor.  The 
intensity,  after  being  divided  by  the  gas  density,  still  shows  a  decrease  with 
increasing  H/C  ratio. 

A  final  expression  was  formulated  by  dividing  the  peak  intensity  by  the 
number  of  carbon-to-carbon  bonds  in  the  fuel  molecule.  For  three  of  the 
fuels,  propane,  ethylene,  and  propylene,  this  gave  a  linear  relationship  when 

plotted  against  pressure.  Figure  35  shows  that  the  value  of  - - - 

#  C— C  bonds 

rises  linearly  with  pressure,  and  the  curves  lie  on  top  of  each  other.  The  slope 
of  these  lines,  12.5,  is  about  l/5  the  slope  from  the  work  of  Claus.  The 
flames  in  Claus’  study  were,  however,  luminous,  and  the  fuel  was  aviation 


kerosene. 
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Emissivities 


A  logical  method  of  analyzing  radiation  data  is  to  calculate  emissivities 
for  the  different  conditions  studied.  In  this  work,  where  the  flame 
temperatures  and  emitted  fluxes  were  measured,  the  emissivities  were 
calculated  as  follows: 


where 


e  =  emissivity 
E  =  emitted  flux  (W /m2) 

a  =  Stefan-Boltzmann  constant  (5.67  *10~8  W/m2-K4) 

T  =  temperature  measured  at  the  plane  of  the  windows 
and  the  subscript  f  denotes  flame  conditions. 

The  emissivities  calculated  from  Eqn  (46)  showed  considerable  scatter.  As 
Figure  36  shows  for  propane,  no  clear  trends  could  be  seen  for  the  effects  of 
equivalence  ratio  or  pressure.  Similar  results  were  observed  for  ethane  and 
propylene. 

Lack  of  a  simple  relationship  describing  flame  emissivity  is  not  surprising, 
especially  in  light  of  the  complex  nature  of  gaseous  radiation.  For  instance, 
the  radiant  energy  emitted  by  C02  and  H20  is  due  to  narrow  bands  in  the 
infrared  portion  of  the  spectrum,  and  therefore  does  not  provide  a  blackbody 
distribution  (0T4).  Nevertheless,  if  the  gas  emissivity  over  each  wavelength 


My 
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band  were  unity,  then  the  overall  emissivity,  e,  would  scale  as  T4.  Since  it 
does  not,  we  can  conclude  that  the  individual  band  emissivities  are  different 
from  unity  and  temperature  dependent.  Considering  that  the  diameter  of  the 
combustion  chamber  is  only  3.81  cm,  this  is  not  surprising,  for  only  when  the 
optical  path  length  becomes  large  do  the  band  emissivities  approach  unity. 

Hottel’s  charts  were  applied  to  the  results  using  the  measured  flame 
temperature  and  combustion  pressure  and  calculated  equilibrium 
concentrations  of  C02  and  H20.  Properties  and  concentrations  were  assumed 
to  be  uniform  across  the  diameter  of  the  combustor.  The  resulting 
(calculated)  emissivities  for  the  non-luminous  gases  show  a  clear  increase 
with  pressure.  Figure  37  shows  this  behavior  for  ethane  and  also  reveals  a 
strong  decrease  in  emissivity  with  increasing  equivalence  ratio.  The  increase 
in  emissivity  with  pressure  is  logical  because  of  collisional  broadening  of  the 
emission-absorption  bandwidths  for  C02  and  H20.  The  decrease  in 
emissivity  with  downstream  location  is  probably  a  manifestation  of  increasing 
temperature  as  '^mbustion  approaches  completion. 

Figure  38  also  shows  emissivities  calculated  using  Hottel’s  charts  and 
indicates  a  decrease  in  emissivity  with  increasing  distance  from  the 
flameholder.  This  behavior,  coupled  with  the  rise  in  temperature  with 
downstream  distance,  explains  why  the  profiles  of  radiation  intensity  show 
little  variation  with  axial  location.  The  variations  in  temperature  and 
emissivity  have  competing  effects,  which  serve  to  flatten  the  flame  radiation 
profiles.  As  Figures  30  through  32  show,  the  flame  temperature  generally 
rises  up  to  x/D  =  4.  The  band  emissivities  are  known  to  fall  with  rising  gas 
temperature  in  this  temperature  range  [Siegel  and  Howell]. 
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Finally,  the  measured  values  for  radiation  intensity  were  compared  to 
values  using  the  calculated  emissivities.  This  allows  for  estimation  of 
radiation  intensity  via  the  equation 


SUMMARY  AND  CONCLUSIONS 


This  concludes  the  study  of  the  effects  of  pressure,  temperature,  and  fuel 
type  on  soot  threshold  equivalence  ratios  and  non-luminous  flame  radiation. 
To  summarize,  the  pressure  ranged  from  0.1  MPa  to  0.8  MPa,  flame 
temperatures  were  as  low  as  1600K  and  as  high  as  2400K,  and  a  variety  of 
paraffins  and  olefins  as  well  as  cyclopropane  were  burned.  All  flames  were 
completely  premixed  and  either  non-luminous  or  just  sooting. 

With  regard  to  soot  thresholds,  this  investigation  has  revealed  several 
important  findings. 

—  Higher  pressures  give  lower  soot  threshold  equivalence  ratios  (<£c).  The 
effect  diminishes  as  pressure  increases. 

—  The  two-step  model  for  soot  precursor  evolution  is  valid  for  wide  ranges  of 
pressures,  flame  temperatures,  and  fuels,  i.e. 


[ XoH  P  ]T 

accurately  describes  incipient  sooting  behavior  over  the  range  of 
conditions  listed  above  and  for  the  fuels  tested  (CH4,  C2Hfl,  C3H8,  C2H4, 
C3H6,  C4H10  and  cyclopropane). 


—  Measurement  of  flame  temperature  at  the  location  of  sooting  is 
imperative;  adiabatic  flame  temperatures  fail  to  explain  the  results. 


—  The  effect  of  fuel  type  on  incipient  soot  formation  appears  negligible. 
This  conclusion  is  based  on  an  analysis  using  the  measured  temperature 
at  that  location  where  soot  is  first  observed  to  form. 

—  That  incipient  soot  behavior  is  independent  of  fuel  type  implies  that  a 
single  chemical  mechanism,  and  single  precursor  species,  is  responsible  for 
the  onset  of  sooting. 

The  second  portion  of  this  investigation  revealed  the  effects  of  equivalence 
ratio,  pressure  and  fuel  type  on  total  non-luminous  flame  radiation. 
Important  features  of  this  portion  of  the  study  were  the  use  of  premixed, 
pure  hydrocarbon  fuels.  As  a  consequence,  the  effects  of  atomization, 
hydrocarbon  blending  and  fuel-air  mixing  were  absent.  Flame  radiation 
could  therefore  be  examined  as  a  separable  function  of  individual  physical 
and  chemical  properties  without  the  complications  associated  with  liquid- 
fuels.  Several  trends  were  apparent. 

—  Radiation  intensity  increased  linearly  with  pressure  up  to  0.8  MPa.  The 
greatest  rate  of  increase  was  for  <f>  —  1.00,  followed  by  4>  =  1.25  and  <j>  = 
0.80. 

—  Radiation  intensity  was  approximately  equal  for  <t>  =  0.80  and  <f>  =  1.25 
and  20%  higher  for  4>  —  1.00. 

—  Radiation  intensity,  normalized  by  the  number  of  carbon-to-carbon  bonds 

in  the  fuel,  — - - ,  scales  with  pressure. 

#  C-C  bonds 

—  Radiation  intensity  decreases  as  the  hydrogen-to-carbon  (H/C)  ratio  of 
the  fuel  increases  under  all  conditions,  which  agrees  with  previous  liquid- 
fueled  studies. 
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N  —  Empirical,  non-luminous  emissivities  increase  with  pressure  and  decrease 

with  axial  location  and  equivalence  ratio. 

—  Emissivities  calculated  from  total  radiation  and  temperature 
measurements  showed  no  trends,  indicating  gas  band  emissivities  of  less 
than  one  for  the  temperatures  and  geometry  involved. 

The  results  of  this  study  suggest  a  follow  up  study  be  conducted  to 
investigate  the  pressure  and  temperature  dependence  of  the  incipient  soot 
behavior  of  mixtures  of  single  component,  vaporized  fuels.  An  investigation 
into  the  d>c  behavior  of  mixtures  would  have  obvious  applications  to  the  new 
class  of  "endothermic"  fuels,  not  only  because  of  the  feasibility  of  modeling 
the  decomposition  products  of  these  fuels  as  a  mixture  of  hydrogen  and 
another  hydrocarbon,  but  also  because  current  plans  call  for  burning 
endothermic  fuels  as  gases  premixed  with  air.  A  study  of  this  type  would 

1  ■  -fc 

•  show  whether  the  mechanism  of  incipient  soot  formation  is  independent  of 

fuel  type.  In  addition,  it  would  have  the  added  advantage  of  providing 
design  engineers  with  values  for  4>c  versus  pressure  and  temperature  so  that 
they  could  design  soot-free  combustors.  Such  combustors  would  have  a 
smaller  radiation  signature,  a  characteristic  of  obvious  importance  to  the  Air 
Force. 
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Appendix  A.  Data  Tables 


Table  Al.  Fuel  Properties. 


Methane 

Ethane 

Propane 

n-Butane 

Ethylene 

Propylene 

1-Butylene 

Cyclopropane 


Formula  M(g/gmol) 


CH4 

C2Hfl 

C3Hg 

C<H10 

C2H4 

c3h6 

C4Ha 

C3Ha 


Hf°(cal/mole) 


-17890 

-20240 

-24820 

-29810 

12500 

4880 

280 

12800 


3 


Table  A 2.  Critical  Equivalence  Ratios. 


Tin(K)  P(MPa) 


Methane 

CH4 


Ethane 

C2H6 


Propane 

C3H8 


rnmmmmmmmmmm 


104 


Fuel 


Tin(K)  P(MPa)  4>c  4 


n-Butane 

c4h10 

Ethylene 

C2H4 


700 

500 

600 

700 


0.2  133 

0.3  1.12 

0.3  1.19 

0.5  109 

0.6  1.07 

0.3  1.19 

0.4  1-15 

0.6  1.06 

0.1  1.37 

0.2  1.55 

0.3  1-37 

0.4  133 

0.5  1-28 

0.6  1-26 

0.7  1-25 

0.8  123 


0.92 

0.78 

0.79 

0.73 

0.71 

0.79 

0.77 

0.71 

1.11 

1.03 

0.91 

0.89 

0.85 

0.84 

0.83 

0.82 


Propylene 

700 

0.1 

0.2 

1.70 

1.72 

1.13 

1.15 

C3HQ 

0.3 

1.48 

0.99 

0.4 

1.31 

0.87 

0.5 

1.29 

0.86 

0.6 

1.27 

0.85 

1-Butylene 

700 

0.1 

0.2 

1.22 

1.12 

0.81 

0.75 

c4h8 

0.3 

1.05 

0.70 

Cyclopropane 

700 

0.1 

0.2 

1.34 

1.18 

0.89 

0.79 

CjHg 

0.3 

1.11 

0.74 

0.4 

1.09 

0.73 

105 


106 


3 


Fuel 

Tit(K) 

P(MPa) 

TrfK) 

10® 

Tf 

Xoh‘107 

n-Butane 

700 

0.2 

1550 

645 

5 

c4h10 

0.3 

2250 

444 

9810 

Ethylene 

500 

0.3 

1580 

633 

8 

C2H4 

0.5 

2300 

435 

11304 

0.6 

2140 

467 

4010 

600 

0.3 

1750 

571 

83 

0.4 

1930 

518 

536 

0.6 

2150 

465 

4670 

700 

0.1 

1370 

730 

- 

0.2 

1580 

633 

4 

0.3 

1610 

621 

8.8 

0.4 

1730 

578 

40 

0.5 

2000 

500 

566 

0.6 

2080 

481 

1130 

0.7 

2200 

455 

2638 

0.8 

2430 

412 

12400 

Propylene 

700 

0.1 

1310 

763 

_ 

c3h6 

0.2 

1460 

685 

1 

0.3 

1965 

509 

335 

0.4 

2050 

488 

987 

0.5 

2065 

484 

1040 

0.6 

2140 

467 

1830 

1-Butylene 

700 

0.1 

1460 

685 

2 

C4H8 

0.2 

1590 

629 

854 

0.3 

1680 

595 

1628 

Cyclopropane 

700 

0.1 

1460 

685 

1 

c3h8 

0.2 

1680 

595 

39 

0.3 

1870 

535 

371 

0.4 

2060 

485 

2110 

|V,T>W«ST' 
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Table  A4.  Adiabatic  Flame  Temperatures  and  Equilibrium  OH  Mole  Fractions. 


Tin(K)  P(MPa)  T*d(K)  —  Xoh'IO7 

1  *d 


108 


Tin(K)  P(MPa)  T^(K)  xOH107 


Table  A5.  (f>c  and  ipc  vs  Measured  Flame  Temperature. 


Fuel  T; 


Methane 

CH. 


Ethane 

C2H8 


Propane 

c3h8 


<Pc 

„  i  11/21 

T  ,  <t>c  T  H  . 

p  ln  ^  '  d  '  n  h 

4 

a 

XOH 

r  J  XOH  P  C 

XOH 

22.6 

23.3 

22.3 

20.8 

21.1 

20.5 

18.1 

20.4 

17.8 

16.1 

18.7 

15.9 

16.0 

17.8 

15.8 

10.5 

20.0 

19.1 

17.7 

18.3 

17.4 

15.5 

16.0 

15.1 

22.3 

22.8 

21.9 

18.5 

19.1 

18.2 

16.4 

16.9 

16.0 

15.8 

16.4 

15.5 

25.7 

26.3 

25.4 

20.6 

21.1 

20.2 

10.0 

20.4 

10.6 

18.1 

18.7 

17.8 

17.3 

17.8 

16.0 

15.8 

16.4 

15.5 

21.6 

22.1 

21.2 

10.0 

19.5 

18.6 

18.4 

18.9 

18.0 

15.0 

16.4 

15.5 

24.3 

24.8 

24.0 

21.7 

22. 

21.4 

10.0 

20.4 

10.6 

18.6 

19.1 

18.3 

16.2 

16.7 

15.0 

15.8 

16.3 

15.5 

23.4 

23.9 

23.0 

19.0 

20.4 

10.6 

18.6 

19.1 

18.2 

17.1 

17.6 

16.7 

23.8 

24.3 

23.5 

20.4 

20.9 

20.1 

18.0 

18.5 

17.6 

16.6 

17.0 

16.2 

14.9 

15.4 

14.6 

i.1 


Table  A6.  <f> .  and  i£c  vs  Adiabatic  Flame  Temperature. 


Or  T  Or  T 

TiJK)  P(MPa)  In  — •  -  In  — — ± 

XoH  ?  XOH  “ 


Methane 

CH. 


Ethane 

C2Ha 


Propane 

c3h8 


700 

0.4 

18.5 

18.2 

0.5 

17.0 

16.7 

0.6 

16.4 

16.1 

0.7 

15.8 

15.5 

0.8 

15.4 

15.1 

500 

0.4 

16.2 

15.9 

0.5 

16.0 

15.7 

0.6 

14.7 

14.3 

600 

0.3 

18.5 

18.1 

0.5 

16.6 

16.2 

0.7 

16.0 

15.7 

0.8 

15.8 

15.5 

700 

0.3 

17.6 

17.3 

0.4 

16.7 

16.3 

0.5 

16.4 

16.1 

0.6 

15.8 

15.4 

0.7 

15.3 

15.0 

0.8 

15.1 

14.7 

400 

0.5 

18.6 

18.2 

0.6 

17.0 

16.6 

0.7 

16.5 

16.1 

0.8 

15.1 

14.8 

500 

0.3 

19.2 

18.8 

0.4 

18.3 

17.9 

0.5 

17.5 

17.1 

0.6 

16.7 

16.4 

0.7 

15.7 

15.4 

0.8 

15.1 

14.8 

600 

0.2 

19.7 

19.4 

0.3 

17.3 

16.9 

0.4 

16.3 

16.0 

0.5 

15.2 

14.9 

0.6 

14.9 

14.5 

0.7 

14.7 

14.4 

0.8 

14.4 

14.0 

700 

0.2 

19.5 

19.1 

0.3 

17.5 

17.1 

0.4 

15.7 

15.3 

0.5 

15.6 

15.2 

0.6 

14.9 

14.6 

0.7 

14.6 

14.3 

0.8 

14.3 

14.0 

8S5 

5? 

Cis-iC 

g 

s 

_  _  _ _ ^ 

1 

III 

Fi 

i  J 

i 

n-Butane 

c4h10 


Ethylene 

C2H< 


Propylene 

C3H6 


1-Butylene  700 

c4h8 


Cyclopropane  700 

C3H8 


w 


VWVV>W?VJirw;<,VVWV,ir.i 


0.80 

P  = 

11.3 

0.4  MPa 

11.5 

11.5 

17.5 

1.00 

10.5 

11.9 

11.9 

18.3 

1.25 

9.25 

11.5 

11.5 

13.7 

0.6  MPa 

0.80 

13.8 

16.8 

15.0 

20.9 

1.00 

15.7 

17.6 

18.1 

25.1 

1.25 

14.0 

16.4 

16.7 

22.7 

0.8MPa 

0.80 

17.1 

18.8 

18.1 

1.00 

20.6 

22.0 

22.0 

- 

1.25 

17.1 

18.8 

18.1 

- 

Temperature  at  Point  of  Maximum  Intensity  (K) 

Ethane  Propane  Ethylene  Propylene 


0.4  MPa 

0.80 

1560 

1670 

_ 

1770 

1.00 

1680 

1870 

- 

1980 

1.25 

1860 

1870 

- 

1960 

0.6  MPa 

0.80 

1590 

2040 

1730 

1.00 

1980 

2150 

- 

2100 

1.25 

2070 

2080 

- 

2180 

0.8MPa 

0.80 

1960 

1970 

_ 

1.00 

2060 

2210 

- 

- 

1.25 

2120 

2180 

- 

- 

V.  A  v  -V  s.  s. 
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Table  A8.  Ethane:  Intensity,  Temperature,  and  Emissivity  Profiles. 
Radiant  Intensities  (kW/m2,sr) 


x/D  =  1 


c/D  —  2 


x/D  =  3  x/D  =  4 


P  =  0.4  MPa 

0  =  0.80 

4>  =  1.00 

0  =  1.25 

11.3 

9.42 

7.68 

9.77 

10.5 

8.38 

7.85 

9.95 

9.07 

7.33 

9.77 

9.25 

P  =  0.6  MPa 

<f>  =  0.80 

13.8 

13.4 

11.3 

9.77 

0  ==  1.00 

14.3 

15.7 

15.0 

14.5 

<f>  =  1.25 

11.9 

13.3 

14.0 

13.9 

P  =  0.8MPa 

0  =  0.80 

16.1 

17.1 

15.0 

13.3 

0  =  1.00 

19.5 

20.6 

20.4 

20.0 

0  =  1.25 

14.7 

15.9 

17.1 

16.8 

Emissivities 


116 


m 

-yi 

••a 


•'•I 


Table  Afl.  Propane:  Intensity,  Temperature,  and  Emissivity  Profiles. 
Radiant  Intensities  (kW/m2,sr) 


x/D  =  1 

x/D  =  2 

x/D  =  3 

x/D 

P  =  0.4  MPa 

<t>  =  0.80 

11.5 

11.5 

10.5 

10 

<f>  =  1.00 

11.9 

11.9 

11.2 

11 

4>  =  1.25 

9.95 

11.5 

11.2 

10 

P  =  0.6  MPa 

<f>  =  0.80 

16.8 

16.1 

14.7 

14 

4>  =  l.oo 

16.1 

17.6 

16.8 

16 

4>  =  1.25 

14.0 

16.1 

16.4 

16 

P  =  0.8MPa 

<t>  =  0.80 

18.1 

18.8 

17.5 

16 

<t>  =  1.00 

18.1 

22.0 

22.0 

21 

<t>  =  1.25 

16.8 

18.8 

18.3 

18 

Temperatures  (K) 

x/D  =  1 

x/D  =  2 

x/D  =  3 

x/D 

P  =  0.4  MPa 

(f)  =  0.80 

4>  =  l.oo 

4>  =  1.25 

i: 

l’ 

li 

P  =  0.6  MPa 

<}>  =  0.80 

21 

4>  =  1.00 

2: 

t >  =  1.25 

li 

P  =  0.8MPa 

<j>  =  0.80 

li 

<t>  =  1.00 

2 

4>  =  1.25 

2 

Table  All.  Propylene:  Intensity,  Temperature,  and  Emissivity  Profiles. 

Radiant  Intensities  (kW/m2-sr) 


=  0.4  MPa 


=  0.6  MPa 


=  0.4  MPa 


=  0.6  MPa 


x/D  =  1 


x/D  =  2 


c/D  =  3  x/D  =  4 


17.5 

14.2 

11.4 

9.77 

14.0 

16.4 

17.8 

18.3 

11.9 

15.4 

16.4 

16.4 

20.9 

18.1 

17.1 

14.5 

23.7 

25.1 

24.4 

23.7 

20.2 

21.6 

22.3 

22.7 

Temperatures  (K) 

T>  =  1 

x/D  =  2 

x/D  =  3 

x/D  = 

1770 

1840 

1840 

1840 

1740 

1920 

1920 

1980 

1670 

1740 

1800 

1960 

1730 

1810 

1960 

2050 

1950 

2100 

2150 

2210 

1900 

2070 

2180 

2180 
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Table  A12.  Empirical  Emissivities  Using  Measured  Flame  Properties. 
Fuel  x/D  =  2  x/D  =  3  x/D  =  4 


Ethane 


P  =  0.4  MPa 


<t>  =  0.80 

.064 

.069 

.061 

<f>  =  1.00 

.074 

.066 

.058 

<t>  =  1.25 

.070 

.058 

.053 

P  =  0.6  MPa 

4>  =  0.80 

.086 

.083 

.070 

<f>  =  1.00 

.077 

.068 

.062 

<t>  =  1.25 

.063 

.060 

.060 

P  =  0.8MPa 

4>  =  0.80 

.083 

.080 

.076 

<t>  =  1.00 

.086 

.076 

.071 

4>  =  1.25 

Propane 

.075 

.072 

.068 

P  =  0.4  MPa 

<t>  =  0.80 

.066 

.066 

.059 

<f>  =  1.00 

.061 

.061 

.058 

4>  =  1.25 

.053 

.053 

.053 

P  =  0.6  MPa 

<t>  =  0.80 

.068 

.068 

.068 

< f)  =  1.00 

.064 

.064 

.064 

4>  =  1.25 

.066 

.063 

.059 

P  =  0.8MPa 

<t>  =  0.80 

.084 

.084 

.075 

4>  =  l.oo 

.074 

.074 

.069 

<t>  =  1.25 

.068 

.068 

.068 
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Appendix  B.  Journal  Publications 

The  following  journal  publications  have  resulted  from  this  work,  have 
been  submitted,  or  are  in  preparation. 


1.  "The  Influence  of  Pressure  and  Temperature  on  Incipient  Soot 
Formation  in  Premixed  Flames,"  W.G.  Cummings,  EH,  P.E.  Sojka,  and 
A.H.  Lefebvre,  submitted  to  Combustion  and  Flame  (July,  1987). 

2.  "Measurements  of  Thermal  Radiation  in  Pressurized  Turbulent 
Flames,"  W.G.  Cummings,  IE,  A.H.  Lefebvre,  and  P.E.  Sojka, 
submitted  to  ASME  Journal  of  Engineering  for  Gas  Turbines  and  Power 
(July,  1987). 
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Appendix  C.  Professional  Personnel 

Three  people  have  been  associated  with  this  contract:  A.H.  Lefebvre, 
Reilly  Professor  of  Combustion,  P.E.  Sojka,  Assistant  Professor  of  Mechanical 
Engineering,  and  W.G.  Cummings,  EH.  Mr.  Cummings’  M.S.M.E.  (Master  of 
Science  in  Mechanical  Engineering)  thesis  was  a  direct  result  of  the  work 
performed  under  this  contract.  The  title  was  "Soot  Thresholds  and  Flame 
Radiation  in  High  Pressure,  Gaseous  Fuel  Flames."  Mr.  Cummings  degree 
was  conferred  December  12,  1886. 


Appe ndix  D.  Inter ac tio ns 


Several  technical  interactions  have  occurred  during  this  contract.  They 
are  divided  into  two  parts.  Critical  soot  equivalence  ratio  and  flame 
radiation. 

Critical  Soot  Equivalence  Ratio 

W.G.  Cummings,  ID,  P.E.  Sojka,  and  A.H.  Lefebvre,  "The  Influence  of 
Pressure  and  Temperature  on  Incipient  Soot  Formation  in  Premixed  Flames," 
Central  States  Section/Combustion  Institute  Spring  Technical  Meeting, 
Argonne  National  Lab,  Chicago  (May  1987). 

W.G.  Cummings,  HI,  P.E.  Sojka  and  A.H.  Lefebvre,  "The  Effects  of 
Temperature,  Pressure  and  Fuel  Type  on  Flame  Radiation  and  Soot 
Formation  in  Gaseous  Hydrocarbon  +  Air  Flames,"  Central  States 
Section/Combustion  Institute  Spring  Technical  Meeting,  NASA  LeRC, 
Cleveland  (May  1986). 

Radiation 

W.G.  Cummings,  m,  P.E.  Sojka,  and  A.H.  Lefebvre,  "Measurements  of 
Thermal  Radiation  in  Pressurized  Turbulent  Flames,"  ASME  1987  National 
Heat  Transfer  Conference,  Pittsburgh  (Aug.  1987). 

W.G.  Cummings,  HI,  P.E.  Sojka,  and  A.H.  Lefebvre,  "Flame  Radiation  in  a 
Gas  Turbine  Combustor,"  SAE  Paper  No.  85208,  SAE  International  Fuels 
and  Lubricant  Meeting,  Tulsa  (Oct.  1985). 

W.G.  Cummings,  HI,  P.E.  Sojka,  and  A.H.  Lefebvre,  "Flame  Radiation  in  a 


